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PERFORMANCE  COMPARISON  OF  TROPOSPHERICPROPAGATION 
MODELS:  RAY-TRACE  ANALYSIS  RESULTS 
USING  WORLDWIDE  TROPOSPHERIC  DATABASES 


1.  INTRODUCTION 

The  physics  of  atmospheric  propagation  in  communication  systems  is  affected  by  ever-changing 
meteorological  conditions  in  the  atmosphere  and  complex  boundary  conditions  on  the  ground.  Whereas 
forecasting  meteorological  conditions  in  small  areas  is  by  itself  difficult,  if  not  impossible,  prediction  of 
propagation  adds  another  dimension  of  difficulty.  This  is  understandable  because  the  boundary  conditions 
of  the  problem  are  complex  and  constantly  changing,  and  the  solution  process  in  itself  is  not  simple 
whether  one  uses  a  wave-optics  or  a  geometrical-optics  approach.  A  great  deal  of  interest  has  been 
focused  recently  in  the  areas  of  signature  analysis,  classification,  and  modeling  of  airbome/spacebome 
synthetic  aperture  radar  (SAR)  data,  topography,  telemetry/command  data,  tracking,  and  accurate 
determination  of  aircraft  or  spacecraft  position  and  orbit  parameters  as  well  as  geophysical  parameters. 
An  important  requirement  within  most  of  those  application  areas  is  the  geometric  and  radiometric 
calibrations  of  airborne  or  spaceborae  data.  In  airborne  and  spacebome  communication  systems,  the 
dynamic  properties  of  the  aircraft  (or  spacecraft)  and  the  atmospheric  turbulence  can  produce  large 
motion  errors,  which  introduce  additional  geometric  and  radiometric  distortions  in  the  data.  Errors  from 
tropospheric  effects  generally  have  been  neglected  during  the  calibration  process  of  those  spacebome  data 
since  other  sources  of  errors  (i.e.,  Doppler  shifts,  ephemeris,  geodetic  systems,  spacecraft  motion,  etc.) 
have  produced  much  larger  errors.  The  performance  of  low-angle  microwave  propagation  over  the  Earth's 
surface  is  the  mode  of  tropospheric  refraction  that  affects  most  types  of  radar  and  navigational  systems 
like  the  microwave  landing  systems  used  at  airports,  as  well  as  line  of  sight  (LOS)  and  mobile  radios. 

Tropospheric  refraction  produces  two  main  effects  on  radio  waves — angular  bending  and  time 
delay.  The  angular  bending  is  due  primarily  to  the  change  of  the  index  of  refraction  with  the  height  of  the 
atmosphere.  The  time  delay  occurs  primarily  because  the  index  of  refraction  is  greater  than  unity,  thus 
slowing  the  speed  of  the  radio  wave,  and  to  a  lesser  extent,  because  of  the  lengthening  of  the  path  by 
angular  bending.  The  correction  measures  of  these  two  tropospheric  propagation  problems  have  been 
proposed  based  on  model-based  or  empirical  measurements  in  terms  of  frequencies  and  climatic  regions 
over  the  years  [1-6]  to  reduce  propagation  errors  in  telecommunication  systems  links.  Many  tropospheric 
models  have  been  proposed  since  the  early  1950s.  Only  five  models  are  selected  in  this  study  (Hopfield, 
Goad,  Blake,  exponential,  and  Cains  with  Case  1)  because  other  models  are  not  directly  related  to  this 
project  or  not  good  enough  to  generate  an  acceptable  range  of  performance.  Little  information  exists  on 
Earth-space  propagation  at  lower  elevation  angles  (<  5°).  However,  important  details  relating  to  the 
properties  of  the  received  signal  are  generally  less  certain.  These  details  might  include  signal  amplitude, 
delay  times  between  different  paths,  and  individual  angles  of  arrival  under  multipath  conditions. 
Meteorological  uncertainties  severely  limit  the  usefulness  of  models  of  existing  microwave  propagation 
passing  through  low  atmosphere  specifically  in  the  presence  of  precipitation  and  in  the  low  elevation 
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angles.  Many  propagation  problems  on  the  LOS  links  arise  from  the  occurrence  of  anomalous  departures 
from  the  normal  value  in  the  vertical  gradient  of  the  refractive.  This  value  itself  will  vary  slowly  with 
season,  time  of  day,  location,  and  the  standard  gradient  in  refractivity  often  being  quoted  as  -40  N- 
unit/km,  corresponding  to  a  4/3  Earth  [6].  The  propagation  effects  that  are  prevalent  when  radio  waves 
traverse  the  atmosphere  manifest  themselves  as  refractive  bending,  time  delays,  Doppler  errors,  rotation 
of  the  plane  of  polarization  (Faraday  effect),  dispersion  effects,  and  attenuation.  The  atmospheric  radio 
refraction  effects  in  the  tropospheric  region  cause  an  extra  time  delay  in  transmission  of  the  signal  and  an 
increase  in  the  elevation  angle  measured  by  the  antenna  system.  In  other  words,  there  are  two  types  of 
errors:  errors  in  measuring  distance  by  means  of  timing  the  transit  of  radio  signals  between  two  points 
(known  as  range  errors)  and  errors  in  estimating  the  elevation  angle  of  a  target  by  means  of  measuring  the 
angle-of-arrival  of  radio  signals  from  the  target  or  spacecraft  (known  as  elevation  angle  errors)  [7].  The 
emphasis  here  is  concentrated  on  the  physical  phenomena  in  the  atmosphere  using  empirical  data  rather 
than  on  building  models  or  analyses  based  on  models  developed. 

New  systems  that  operate  at  low  elevation  angles  require  improved  accuracy  in  range  errors  and 
angle-of-arrival  errors.  An  approach  for  obtaining  more  accurate  angle  and  range  error  corrections  is  to 
use  calibration  sources  such  as  the  limb  of  the  Sun,  radio  sources,  or  satellites  for  angle  error  cases. 
Neither  radio  sources  nor  satellites  are  suitable  for  range  error  calibration  since  radio  sources  are 
collected  passively,  and  the  true  range  of  a  satellite  is  not  generally  known  to  a  sufficient  degree  of 
accuracy  to  be  of  value.  In  the  formulation  of  these  problems,  an  idealized  model  of  a  time  independent, 
spherically  stratified  nonionized  atmosphere  with  an  index  of  refraction  that  monotonically  decreases 
with  increasing  altitude  is  adopted  for  a  newly  proposed  model.  The  new  tropospheric  model  is  proposed 
here  with  five  other  potential  tropospheric  propagation  models  [1-6]  developed  over  the  last  3  decades; 
they  are  presented  for  comparison  purposes  to  correct  lower  angle  propagation  errors  such  as  time  delays 
or  range  and  angle  errors.  Important  details  relating  to  the  properties  of  the  received  signals  are  generally 
less  certain;  such  details  might  be  signal  amplitude,  delay  times  between  different  paths,  and  individual 
angles  of  arrival  under  multipath  conditions. 

Tbere  are  many  variables  that  influence  Earth-space  propagation.  They  fit  into  the  broad  categories 
of  frequency,  space,  and  time.  Propagation  effects  depend  on  all  these  variables  and  are  quantified 
through  measurements  or  modeling.  In  some  cases,  measured  values  can  be  applied  directly.  In  other 
words,  data  might  be  available  for  the  particular  frequency,  elevation  angle,  and  climate  zone  for  a 
proposed  system.  More  commonly,  measurements  are  used  in  combination  with  theoretical  calculations 
from  fundamental  propagation  physics  to  develop  models  that  explain  the  variations  and  that  can  be 
evaluated  for  specific  situations.  In  fact,  models  usually  are  used  to  predict  average  behavior  while 
measurements  over  some  periods  reveal  the  year-to-year,  season-to-season,  or  day-to-day  variability.  In 
order  to  be  comprehensive,  this  report  covers  the  analysis  results  for  both  selected  leading  models  and 
three  measured  worldwide  climatological  and  diurnal  meteorological  databases.  A  brief  overview  of 
database  characteristics  that  will  be  applied  in  this  study  is  presented  in  Section  2.  Each  leading  model 
with  a  new  proposal  is  introduced  in  Section  3  with  details  of  mathematical  and  physical  principles. 
Results  of  this  study  are  described  in  Section  4.  Conclusions  and  recommendations  are  in  Section  5, 
followed  by  references  and  acknowledgments.  Finally,  figures  and  tables  of  analysis  results  are  included 
in  the  Appendixes. 

2.  DATABASE  CHARACTERISTICS 

Databases  make  use  of  raw  information  provided  by  the  Air  Force,  Navy,  and  National  Oceanic  and 
Atmospheric  Administration  (NOAA).  Software  developed  on  the  project  takes  these  raw  data  and 
formats  them  into  a  standard  Naval  Research  Laboratory  (NRL)  format.  The  software  makes  use  of 
tropospheric  models  that  permit  the  addition  of  tropospheric  refractivity,  grid  number,  and  related 
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statistics  with  height  information.  The  wind  speed,  wind  direction,  and  precipitation  can  be  added  upon 
request  if  the  customer  requires  this  additional  information.  The  resulting  database  of  refractivity  and 
related  statistics  is  created  using  either  EMPRESS  database  software  or  FORTRAN  and  C  on  Sun 
SparcStations.  All  refractivity  and  related  statistical  data  are  stored  in  a  readable  text  ASCII  format  and 
are  also  available  in  Tar  and  VAX  formats.  These  17  variable  outputs  are  available  in  hard  copy  or  as  soft 
copy  on  8-mm  magnetic  tapes. 

2.1.  Data  Sources 

This  task  has  produced  six  databases  using  data  from  the  following  three  government  agencies;  the 
Air  Force,  the  Navy,  and  NOAA.  Results  obtained  with  these  data  will  be  available  to  government 
agencies  and  laboratories  for  further  research  and  modifications.  The  six  data  sources  currently  supported 
are  the  following: 

1.  European  Center  for  Medium-Range  Weather  Forecast  (ECMWF)  from  the  U.S.  Air  Force 

Environmental  Technical  Applications  Center  (ETAC)  at  Scott  Air  Force  Base,  Illinois. 

•  includes  data  averaged  monthly  over  the  years  1981  to  1991  with  2.5°  x  2.5°  grid; 

•  17-layered  data  by  geopotential  height  and  pressure  levels  from  10  to  1000  mbar; 

•  other  data  elements  include  latitude,  longitude,  temperature,  dew  point,  air  density,  and  number 
of  observations  used  to  obtain  mean  and  standard  deviation  values. 

2.  Asheville  Marine  data  from  the  National  Climatology  Data  Center  (NCDC),  NOAA. 

•  includes  upper-air  diurnal  data  over  the  period  of  January  1980  to  June  1993  with  2.5°  x  2.5° 
grid; 

•  40-layered  data  by  geopotential  height  and  predefined  pressure  from  10  to  1000  mbar  over  the 
coastal  lines  and  oceans  around  the  world; 

•  other  data  elements  include  temperature,  dew  point,  height,  number  of  levels,  latitude,  longitude, 
date,  and  time. 

3.  Fleet  Numerical  Meteorological  and  Oceanographic  Center  (FNMOC)  data  from  the  Naval 

Meteorology  and  Oceanography  Command  in  Monterey,  California. 

•  includes  mean  sea-level  and  upper  air  data  every  12  h  from  January  1994  to  February  1996 
with  2.5°  X  2.5°  grid  around  the  world; 

•  17-layered  data  by  geopotential  height  and  predefined  pressure  level  from  200  to  1000  mbar; 

•  other  data  elements  include  air  temperature,  dew  point,  wind  vector  and  speed,  latitude,  longitude 
date,  and  time  with  2.5°  x  2.5°  grid. 

4.  High-Resolution  Analysis  System  (HIRAS)  data  from  ETAC,  Scott  Air  Force  Base,  Illinois 

•  includes  the  monthly  and  6  hourly  averaged  climatology  data  over  the  period  of  July  1988  to 
June  30,  1994  with  2.5°  x  2.5°  grid  around  the  world; 

•  17-layered  data  by  geopotential  height  and  predefined  pressure  level  from  10  to  1000  mbar; 

•  other  data  elements  include  temperature,  relative  humidity,  height,  latitude,  longitude,  date,  and 
time. 

5.  Medium-Range  Forecast  (MRF)  data  from  NCDC,  NOAA 

•  includes  6  hourly  diurnal  meteorological  data  over  the  period  of  January  1,  1991,  to  Dec¬ 
ember  31,  1995,  with  2.5°  x  2.5°  grid  around  the  world; 

•  13-layered  data  by  geopotential  height  and  predefined  pressure  level  from  50  mbar  (24  km 
from  the  surface  in  the  air)  to  Earth  surface; 
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•  other  data  elements  include  air  temperature,  relative  humidity,  wind  vector  and  speed,  latitude, 
longitude,  date,  and  time 

6.  Final  Analysis  (FNL)  Data  for  MRF 

•  includes  6  hourly  diurnal  meteorological  data  over  the  period  of  January  1,  1997,  to  present  with 
1.0°  X  1.0°  grid  around  the  world; 

•  14-layered  data  by  geopotential  height  and  predefined  pressure  level  from  20  mbar  (26.6  km 
from  the  surface  in  the  air)  to  Earth’s  surface; 

•  other  data  elements  include  temperature,  relative  humidity,  total  cloud  cover,  wind  vector  and 
speed,  geopotential  height  and  pressure  vertical  velocity. 

2.2  Accessing  Databases 

Details  of  accessibility  for  each  database  are  referred  to  in  Ref.  8.  HIRAS  and  MRF  databases  are  not 
included  in  Ref.  1,  and  their  configuration  management  is  similar  to  others  and  will  be  provided  upon 
request. 

3.  MODELING 

A  discussion  of  tropospheric  effects  on  radio  waves  can  be  divided  into  two  parts — the  refractivity 
model  and  range  or  angle-of-arrival-error  model.  Many  models  involving  refractivity  or  range  and 
bending  errors  have  been  proposed  during  the  last  several  decades.  It  is  impossible  to  cover  here  all  the 
models  published  thus  far.  Rather,  this  report  concentrates  on  a  few  leading  models  on  both  refractivity 
and  range  or  angle  errors  in  order  to  present  a  feasible  approach  that  is  useful  for  system  implementation 
in  real-world  applications,  as  pointed  out  in  Section  1. 

3.1  Refractivity  Models 

Since  the  introduction  of  refractivity  N  by  Smith  and  Weintraub  in  1953  [9]  as 

N=  (11. 6/Tip  +  4,810  *  e/T\ ,  (1) 

(where  T  is  the  temperature  in  Kelvin,  P  the  pressure  in  mbar,  and  e  the  water  vapor  pressure  in  mbar), 
scientists  and  engineers  have  proposed  numerous  refractivity  profiles  to  understand  the  propagation  path 
in  the  atmosphere  for  the  compensation  of  range  and  angle  error  at  low  elevation  angles.  Angular  bending 
is  due  primarily  to  the  change  in  the  index  of  refraction  with  the  height  of  the  atmosphere.  Time  delay 
occurs  primarily  because  the  index  of  refraction  is  greater  than  unity  in  the  tropospheric  region  to  the 
height  of  about  30  km  from  the  ground,  thus  slowing  the  radio  wave,  and  to  a  lesser  extent,  because  of  the 
lengthening  of  propagation  path  by  angular  bending.  In  general,  refractive  errors  increase  with  decreasing 
elevation  angle  for  a  standard  atmosphere. 

3.1.1.  Effective  Earth  Radius  Model 

The  Earth  radius  model  was  formulated  by  Schelleng  et  al  [10]  and  was  first  used  for  the  LOS 
communications  problems.  It  was  shown  that,  by  assuming  that  the  Earth  has  a  radius  of  about  4/3  that  of 
the  actual  Earth,  radio  wave  rays  could  be  drawn  as  straight  lines.  It  is  evident  that  the  4/3  Earth 
atmosphere  has  about  the  correct  slope  in  the  first  kilometer  or  two  above  the  Earth’s  surface  but 
decreases  rapidly  above  that  height.  From  an  examination  of  many  years  of  /^-profile  data  for  various 
climates,  the  observed  refractivity  distribution  is  more  nearly  an  exponential  function  of  height  than  a 
linear  function  of  height  as  assumed  by  the  4/3  Earth  atmosphere.  One  might  expect  that  refractivity 
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decreases  exponentially  with  height  since  the  first  term  of  Eq.  (1)  involving  pTT  comprises  at  least  70%  of 
the  total  and  is  proportional  to  air  density,  a  well-known  exponential  function  of  height.  It  appears  that 
this  success  is  due  to  the  4/3  Earth’s  model  essentially  being  in  agreement  with  the  average  N  structure 
near  the  Earth’s  surface  that  largely  controls  the  refraction  of  radio  wave  rays  at  small  values  of  elevation 
angle  common  in  tropospheric  communication  systems.  Based  on  numerous  studies,  refractivity  N  may 
be  represented  by  an  exponential  function  of  height  of  the  form: 

N{h)=N^txp{-bh},  (2) 

where  N  is  the  surface  refractivity,  b  the  constant,  h  the  altitude  from  the  surface  in  the  altitude  range  of  1 
to  9  km  above  the  sea  surface  level.  The  effective  Earth  radius  model  works  well  for  propagation  paths  at 
low  altitudes  where  ray  paths  are  within  about  2  km  of  the  Earth’s  surface  but  not  for  those  at  higher 
altitudes.  Further  details  of  Eq.  (2)  can  be  broken  down  into  the  region  as 

N{h)  =  A(s  +  {h-h^)  AN  for  h^<h<h^  +  l,  (3) 

where  -  AN  =  7.32  exp  {0.00557  N^)  and  is  the  surface  height.  (4) 

Equations  (3)  and  (4)  are  based  on  the  effective  Earth’s  radius  concept  in  the  first  kilometer  from  the 
surface.  In  this  atmosphere,  N  is  assumed  to  decay  linearly  with  height  from  the  surface  to  1  km  above 
the  surface  h^  +  l. 


N(h)  =  Ni  exp  {-c  [h-h^- 1)},  for  +  1  <  /z  <  9  km,  (5) 

where  c  =  {l  /(8  -  /i ln(  / 1 05) ,  (6) 

and  A/j  is  the  value  of  at  1  km  above  the  surface.  Above  the  altitude  of  9  km,  where  less  than  10%  of 
the  total  bending  occurs,  a  single  exponential  decrease  of  N  may  be  assumed.  The  coefficients  in  the 
exponential  expression: 


A^(/i)  =  105exp  {-0.1424(/i-9)},  for  A  >  9km  (7) 

were  derived  by  the  Rocket  Panel  data  [11]. 

The  three-part  model  of  the  atmosphere  expressed  by  Eqs.  (3),  (5),  and  (7)  has  the  advantage  of  the 
effective  Earth’s  radius  model  approach,  particularly  for  such  applications  as  point-to-point  radio  relaying 
over  distances  up  to  100  miles  where  the  radio  energy  is  generally  confined  to  the  first  kilometer  and 
being  in  reasonably  good  agreement  with  the  average  //-structure  of  the  atmosphere.  Note  that  the  4/3 
Earth  model  with  its  constant  decay  of  39.2  N  units  per  kilometer  would  be  a  poor  representation  of  the 
maximum  profile,  which  decreases  over  66  //-units  in  the  first  kilometer.  This  implies  that  the  4/3  Earth 
model  closely  represents  the  slope  of  the  minimal  N^  profile  over  the  first  kilometer  but  then  decreases 
too  rapidly  with  height. 

3.1.2.  Exponential  Model  [1, 12] 

If  it  is  supposed  that  the  Earth’s  atmosphere  consists  entirely  of  isothermal  ideal  gas,  the  tropospheric 
refractivity  profile  is,  as  is  well  known,  calculated  from  the  state  of  the  ideal  gas  equation.  If  the 
refractivity  N^  at  the  height  is  equal  to 
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N(h)  =  N,  exp  {-  {gM /  RT)  {h  -h,)},  (8) 

or  simply 

/V(/i)=A^,exp{-Q(/i-/i,)},  (9) 

where  g  is  the  acceleration  of  gravity  (9.80  m/s  2),  M  the  gram  molecular  weight  of  the  air  (29.0  g  ),  /?  the 
gas  constant  (8.3144  J/mol-K),  T  the  absolute  temperature  in  Kelvin,  and 

Cg  =  In  [/Vj  /Ml .0 km) ]  or  =  In  /(M  +  AN)]  ■  (10) 

These  models  of  atmospheric  refractivity  (Eqs.  (8)  and  (9))  are  a  close  representation  of  the  average 
refractivity  structure  within  the  first  3  km.  Further,  the  single  exponential  model  has  the  advantage  of 
being  an  entire  function  and,  therefore,  is  easily  used  in  theoretical  studies.  The  exponential  reference 
atmosphere  is  in  good  agreement  with  the  initial  N  distribution  but  tends  to  give  systematically  low 
values  above  ~3  km.  Therefore,  the  exponential  reference  atmosphere  does  not  appear  to  be  as  good  a 
representation  of  the  two  observed  profiles  as  the  reference  atmosphere,  particularly  above  approximately 
5  km. 

3.1.3.  Hopfield  Model  [2] 

The  actual  atmosphere  is  not  isothermal  nor  is  its  composition  an  ideal  gas,  and  the  state  of  water 
vapor  in  the  atmosphere  is  different  from  the  state  of  an  ideal  gas.  Therefore,  it  is  more  practical  to 
express  the  refractivity  in  the  quartic  term  induced  by  the  dry  air  or  the  water  vapor  separately  as  follows: 

N{h)  =  Nj  +  N^,  (11) 

where  Nj  is  the  dry  refractivity  and  the  wet  refractivity  represented  by 


Nd  =  kd(hoa-hy,  (12) 

N^=k^(ho^-hy,  (13) 

with  =  l/[/io4  “  kyv-  ^Qd  height  of  the  order  of  40  km,  and  the  wet 

height  of  the  order  of  12  km.  Note  that  if  the  refractivity  as  a  function  of  height  is  represented  by  an 
exponential,  it  is  not  integrable  in  closed  form,  where  if  it  has  the  form  as  in  Eqs.  (11)  through  (13),  it  is 
integrable.  The  representation  of  the  dry  and  wet  terms  of  refractivity  by  quartic  Eqs.  (11)  through  (13) 
gave  good  agreement  with  range  error  and  Doppler  data  above  6°  elevation  angles.  Note  that  if  monthly 
or  weekly  averages  of  the  refractivity  are  used  and  the  accuracy  is  not  so  important,  this  model  is  useful 
and  practical. 

3.1.4.  Modified  Exponential  Model 

It  has  been  seen  that  the  observed  refractivity  distribution  is  more  nearly  an  exponential  function  of 
height  than  a  linear  function,  as  assumed  by  the  effective  Earth’s  model.  The  exponential  decrease  of  the 
refractivity  N  with  height  is  sufficiently  regular  as  to  permit  a  first  approximation  of  average  refractivity 
N  structure  from  surface  condition  alone.  Consider  that 


N{h)  =  N^exp(-h/H), 


(14) 
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where  H  is  a  scale  (or  reference)  height  appropriate  to  the  value  of  N  at  zero  height  N^.  Considering  a 
scale  height  here,  it  is  simply  the  height  at  which  the  value  of  N  (h)  is  1/e  of  under  the  assumption  of 
Eq.  (14),  at  which  the  height  h  is  equal  to  the  scale  height  H.  The  wet  refractivity  is  below  1.0  A^-unit 
in  comparison  with  the  dry  refractivity  with  100  //-units  in  the  neighborhood  of  the  selected  value  of  a 
reference  height  H.  The  ratio  of  dry  to  wet  refractivity  is  approximately  100  at  the  reference  height  H 
where  the  height  from  the  surface  is  one-third  of  the  total  tropospheric  region.  The  refractive  phenomena 
of  bending  and  time  delay  beyond  this  layer  (reference  height)  will  be  limited  since  temperature  and 
humidity  do  not  change  drastically  to  affect  refractive  bending  in  those  extended  areas  such  as 
tropopause,  stratosphere,  stratopause,  free  space,  and  ionosphere  above  1  GHz.  This  coincides  with  the 
fact  that  most  bending  and  refractive  phenomena  occurs  within  this  region  (beneath  the  reference  height) 
from  the  surface  of  the  Earth.  This  implies  that  the  tropospheric  effects  on  the  ray  bending  can  be 
approximated  with  the  reference  height  without  significant  loss  of  any  physical  or  atmospheric  theory. 

As  it  is  well  known,  the  atmospheric  pressure  tends  to  decrease  exponentially  in  accordance  with 
[13, 14] 


P  =  PQQxp(-h/H),  (15) 

where  h  is  the  height  above  a  reference  level  where  the  pressure  is  Pq.  It  is  noted  that  the  scale  height  H, 
however,  is  not  a  constant  as  it  is  a  function  of  temperature  T,  the  average  mass  M  of  the  molecules 
present,  and  the  acceleration  of  gravity  g  as  shown  in  the  Eq.  (8)  by 

H^kTIMg,  (16) 

where  k  is  Boltzmann’s  constant.  The  rate  of  change  of  temperature  with  altitude  in  a  dry  atmosphere  in 
an  adiabatic  state  involving  no  input  or  loss  of  heat  energy  is  given  by 

dT/dh  =  -9.8°  C!/km  (17) 

If  the  actual  lapse  rate  of  the  atmosphere  (rate  of  decrease  of  temperature  with  altitude)  is  9.8°  C/km,  a 
parcel  of  air  that  is  originally  in  equilibrium  with  its  surroundings  and  which  is  then  moved  upwards  or 
downwards  will  tend  to  remain  in  equilibrium  at  the  same  temperature  as  its  surroundings.  Then  the 
parcel  of  air  will  not  be  subject  to  any  restraining  or  accelerating  force.  Such  a  lapse  rate  of  temperature  is 
referred  to  as  neutral.  If  the  actual  lapse  rate  of  atmosphere  is  greater  than  9.8°  C/km,  a  rising  parcel  of  air 
will  tend  to  cool  only  at  the  adiabatic  rate  and  is  warmer  than  its  surroundings.  As  a  result,  it  will  be 
lighter  than  the  air  around  it  and  will  be  accelerated  still  further  upwards. 

In  an  inversion  layer,  temperature  increases  with  altitude,  and  such  a  layer  is  highly  stable.  All 
vertical  motions  are  strongly  inhibited  in  an  inversion  layer,  and  pollution  emitted  below  the  layer  tends 
to  be  confined  below  it.  Also,  if  a  source  of  water  vapor  exists  below  an  inversion  layer,  it  tends  to  be 
confined  below  the  layer,  with  the  result  that  large  decreases  in  index  of  refraction  may  be  encountered  in 
the  upward  passage  through  an  inversion  layer.  Thus,  the  occurrence  of  inversion  layers  has  an  important 
effect  on  low-elevation  angle  Earth-space  communication  paths.  The  decrease  or  change  of  the  water 
vapor  pressure  e  with  height  is  generally  variable  but  may  be  approximately  exponential.  Note  also  that 
the  delay  caused  by  water  vapor  is  considerably  smaller  than  that  for  dry  air  above  3  to  5  km  from  the 
surface,  but  total  water  vapor  content  along  a  path  is  variable  and  not  predictable  with  high  accuracy  from 
the  surface  water  vapor  pressure  or  density.  Therefore  water  vapor  is  responsible  for  a  larger  error  or 
uncertainty  in  the  range  than  in  dry  air  at  lower  atmosphere. 
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This  kind  of  exponential  model  is  widely  applicable  and  is  dependable  when  reliable  climatological 
or  meteorological  data  on  actual  refractivity  profiles  are  applied.  In  this  report,  a  worldwide  2.5°  x  2.5° 
grid  accuracy  is  used.  This  model  provides  the  accuracy  of  less  than  1%  of  root-mean  square  (rms)  error 
from  the  climatology  or  meteorology  data  in  comparison  with  the  accuracy  of  20%  to  30%  of  rms  errors 
for  the  Hopfield  and  other  models.  Therefore,  this  model  approach  has  been  chosen  here  as  the  most 
reliable  and  accurate  in  comparison  with  other  models  for  various  conditions.  The  comparison  and  tested 
results  are  presented  in  Section  4  for  both  spatially,  temporally,  and  geographically  diverse  environmental 
conditions. 

3.1.5.  Complex  Refractivity  Model  [15, 16, 17] 

With  the  current  high  interest  in  millimeter  and  submillimeter  waves,  there  is  a  need  for  a  reliable 
model  to  predict  average  loss  and  delay  effects  from  easily  obtained  meteorological  data.  Such  a  model 
would  find  practical  application  through  conversion  of  basic  climatological  variables  (i.e.,  temperature  T, 
barometric  pressure  P,  relative  humidity  Q)  into  transfer  characteristics  of  a  radio  path.  In  atmospheric 
turbulence,  the  fluctuations  in  T,  P,  and  Q  cause  fluctuations  in  both  the  real  and  imaginary  parts  of  the 
refractive  index.  Such  fluctuations  cause  random  refraction  and  absorption  of  electromagnetic  waves 
passing  through  the  medium.  At  visible  and  radio  frequencies,  the  refractive  index  is  a  relatively  simple 
function  of  T,  P,  and  Q,  and  it  is  fairly  easy  to  express  the  fluctuations  of  the  refractive  index  in  terms  of 
the  fluctuations  of  T,  P,  and  Q,  as  shown  in  Eq.  (1).  However,  for  electromagnetic  radiation  at  millimeter 
and  submillimeter  waves,  the  presence  of  absorption  resonance  causes  both  real  and  imaginary  parts  of 
the  refractive  index  to  depend  on  T,  P,  Q,  and  frequency  in  a  more  complicated  manner.  The  complex 
refractivity  N,  expressed  in  terms  of  measurable  quantities,  provides  that  role.  For  air,  N  consists  of  three 
components: 


N=NQ  +  D(j)  +  jir'(f),  (18) 

namely,  frequency  independent  refractivity  Nq  plus  various  spectra  of  refractive  dispersion  D(f)  and 
absorption  A/"  (f).  The  imaginary  part  of  Eq.  (18)  is  usually  expressed  as  the  specific  power  attenuation  a, 


and  the  real  part  determines  the  phase  delay  P  (with  reference  to  vacuum).  That  is 

a  =  0.1820  JN"  (f)  dB/km,  (19a) 

P  =  0.02096  +  D)  rad/km.  (19b) 

Accordingly,  the  propagation  constant  F  and  the  excess  propagation  delay  time  t  are 

r  =  -4).1151  a  +  X2.096  *  10^/+  p)  1/km,  (20a) 

t  =  ( p/27r^*103  =  3.336  {Nq+  D)  ps/km.  (20b) 


where  ps  denotes  picosecond.  Note  also  that  water  vapor  refractivity  is  about  16  times  more  effective  on  a 
per-molecule  basis  than  dry  air  in  generating  propagation  phenomena  such  as  time  delay,  ray  bending, 
ducting,  scintillation,  etc.  The  absorption  and  dispersion  spectra  are  formulated  from  the  contributions  of 
a  continuum  and  a  liquid  water  extinction  N"  i.e., 

hr'{f)=  l(5F")i+Ar^+Arv 

/ 


ppm 


(21) 
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and 


D(f)=  ppm,  (22) 

/ 

where  SF"  is  the  line  spectra  of  the  absorption  and  SF'  the  refractive  dispersion  with  strength  S  in  units  of 
kilohertz  and  shape  factors  F'  and  F"  in  units  of  (GHz)“k  Since  both  Eqs.  (21)  and  (22)  require  more 
elaboration,  one  can  refer  to  Ref.  15  for  further  details  of  derivations  and  characteristics. 

3.2.  Range  and  Angle  Error  Models 

Many  models  have  been  proposed  to  investigate  the  causes  of  the  observed  errors  in  range  and  angle- 
of-arrival  errors  and  to  determine  the  orbital  elements  of  an  artificial  satellite.  These  models  include  the 
azimuth-elevation  angle  method,  the  radio  interferometer  method,  and  the  laser  ranging  method.  The 
radio  wave  propagation  between  a  ground  station  and  a  spacecraft  is  subject  to  the  bending  of  the 
propagation  path  and  the  decrease  in  the  propagation  velocity  in  the  Earth’s  atmosphere.  These  effects 
cause  systematic  and  random  errors  in  the  range,  range  rate,  and  the  angle-of-arrival  measurements.  The 
errors  caused  in  the  ionosphere  can  be  reduced  to  as  small  as  desired  by  the  use  of  shorter  waves  such  as 
millimeter  or  centimeter  waves,  and  the  errors  in  the  troposphere  are  also  reduced  to  about  an  one-half  or 
two-thirds  if  the  proposed  approach  is  adopted  using  an  available  climatology  database. 

There  are  different  approaches  used  to  compute  range  and  range  errors  that  consider  refraction  effects 
to  improve  range  measurements  by  removing  systematic  bias.  The  range  error  is  composed  of  three  parts: 
the  difference  between  the  curved  length  of  the  propagation  ray  path  R  and  the  true  slant  range  Rq,  mainly 
due  to  the  increase  in  time  necessary  to  travel  over  the  curved  path  R ;  timing  errors  in  the  detection 
system;  and  the  discrepancy  caused  by  the  lowered  velocity  of  propagation  in  a  refractive  medium. 
Propagation  times  for  waves  traveling  between  the  ground  station  and  a  spacecraft  are  longer  than  the 
figures  calculated  for  open  space  for  two  reasons: 

1 .  The  path  does  not  follow  a  straight  line.  The  consequent  increase  in  path  length  is  small  and  can  be 
neglected  except  for  angles  of  elevation  below  5°  to  10°. 

2.  The  radio  wave  velocity  is  slightly  lower  than  it  would  be  in  a  vacuum,  producing  an  apparent 
increase  in  the  length  of  the  path  given  by  the  relation: 

AR  =  J^(n  -  1)  ds,  (23) 

where  s  is  the  curved  abscissa  on  the  path,  and  R  the  distance  of  the  spacecraft,  which  can  be  treated  as 
infinite  for  the  purpose  of  these  calculations  as  atmospheric  effects  only  influence  the  first  few 
kilometers.  As  the  real  path  does  not  deviate  much  from  a  straight  line  with  the  angle  of  elevation  0,  as 
long  as  0  is  greater  than  a  few  degrees,  the  range  error  can  be  approximated  as: 

A/?  =  /o[(n-l)/sin0o]d/i,  (24) 

where  h  is  the  vertical  altitude.  Different  range  error  models  have  been  derived  by  solving  Eqs.  (23)  and 
(24)  in  direct,  linear  approximation  and  discretized  (or  stratified)  approaches.  Note  here  that  no  effort  has 
been  made  to  examine  or  derive  analytical  expressions  of  models.  Rather,  it  introduces  main  model 
equations  to  compare  each  approach.  If  one  wants  to  understand  further  details  of  model  expressions  and 
their  reasonings,  please  refer  to  the  references  provided  for  each  model. 
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3.2.1.  Hopfield  Range  Error  Model  [2] 

Hopfield  assumed  that  the  lapse  rate  of  temperature  with  respect  to  height  is  constant  in  the 
troposphere.  A  lapse  rate  of  6.8°  C/km  was  assumed.  Therefore,  the  functional  form  of  the  dry  refractivity 
in  the  troposphere  becomes  a  quartic.  Hopfield  used  a  linear  approximation  technique  to  solve  Eq.  (24)  as 

Aptro=  lApi,  (25) 

/ 

where  7=1  denotes  dry  component,  1  =  2  wet  component  and 

Ap,=  10-*W„  {-I,+(4/*y(l/3)r?lf-  (2/15)1?  -  (3/4)75. ri„,li(lf  +  (l/2)li) 

+  rL  ll-  (I'3)'i5»l3i  +  P'l5)l|i 

(I3, +(1/2)13,12)  -  dil3i(l3,  -  (l/2)>iL) 

+  (l'2)'i™.  l2((3/4)li+r,?„)  ta[(>5-  +  lT)/(<i„i+l3i)]|,  (26) 

f  1  2  ^  ^ 

and  1  j  =  rj-sin  6, 12  =  rj  cos  6, 13,-  =  - 12  I  ,  r^and  are  distances  from  the  center  of  the  Earth 

to  the  ground  tracking  station  and  to  the  top  layer  of  the  troposphere  (dry  or  wet  component), 
respectively.  Hopfield  adopted  the  dry  height  to  40  km  and  the  wet  height  to  12  km.  Details  of  the 
derivation  should  be  referred  to  Hopfield  [2]. 

3.2.2.  Stratified  Layer  Model  [4, 16, 18] 

If  we  limit  ourselves  with  refraction  in  the  vicinity  of  polar  and  equatorial  regions  where  the  effect  of 
the  Earth’s  magnetic  field  is  an  important  consideration,  it  is  a  convenient  and  valid  approximation  to 
consider  the  atmosphere  as  consisting  of  several  spherically  stratified  layers  within  a  small  segment  of 
Earth  surface  like  a  2.5°  x  2.5°  grid  area  in  the  globe.  If  this  medium  is  slowly  varying  with  height,  it  is 
then  possible  to  assign  numbers  representing  the  mean  value  of  the  refractive  index  at  any  given  time 
(i.e.,  0000,  0600,  1200,  1800  h)  for  each  of  these  layers.  The  bending  of  a  ray  as  it  traverses  each 
successive  layer  is  computed  from  a  form  of  Snell’s  law  that  applies  to  spherically  refracting  surfaces. 
Snell’s  law  for  the  refraction  of  electromagnetic  waves  at  a  plane  interface  between  two  mediums  of 
index  njand  n2  is  given  by 

nj  sin  j I  =  n2 cos  ai ,  (27) 

where  Jj  is  the  angle  of  incidence,  and  ^2  is  the  refraction  angle  in  mediums  nj  and  n2>  respectively. 
Consider  the  small  vicinity  around  the  point  where  the  ray  intersects  the  spherical  boundary  between  the 
two  different  media,  as  Fig.  1  shows.  If  this  region  is  chosen  small  enough  to  be  physically  and 
mathematically  acceptable,  it  may  be  considered  a  plane  and  Snell’s  law  be  assumed  to  apply.  Also,  from 
the  law  of  sines. 


sin  / j/rj  =  sin  (90°  +  cc\)lr2  =  cos  ct j/r2. 


(28) 
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Center  of  Earth 

Fig.  1  —  Progressive  ray  bending  traversing  spherical  atmospheric  layer  stratification 


Also, 


sin  i]  =  (r|/r2)  *  cos  oj 


(29) 


Substituting  Eqs.  (28)  and  (29)  into  Eq.  (27)  yields  Bouger’s  rule, 


n  I  r^  cos  a  j  =  n2r2  cos  02 .  (30) 

The  basic  assumption  that  the  proposed  mathematical  approach  embodies  is  the  following;  the 
atmosphere  is  considered  to  be  stratified  into  m  spherical  layers  of  thickness  and  constant  refractive 
index  n^.  This  type  of  stratification  is  seen  in  Fig.  1  where  Oq  is  the  apparent  elevation  angle,  and  Oq^  is 
the  trae  elevation  angle.  The  general  expressions  for  and  are  given  by 

“m  =  cos-l  /n„  r„) cos  a^_i] ,  (31) 

and 


(32) 
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where  the  radial  distance  is  merely  the  summation  of  the  various  layers  expressed  by 


fm+l  =’-o+^hj. 
j=0 


(33) 


Applying  the  law  of  sines  for  the  direct  path,  it  follows  that 


L  /^0m)sin 

1 

U"®  JJ 

(34) 


where 


=  fQ  +  r^+\  -  2ro  ror^+i  cos 


I  e, 

7=0 


(35) 


and 


0  j  =  7t/2  -  aj  -  ij . 


(36) 


The  refraction  angle  error  Aa^,  which  is  the  difference  between  the  apparent  elevation  angle  and  the  true 
elevation  angle,  can  then  be  determined  from 


(37) 

Similarly,  the  range  error  AR,  which  results  from  the  velocity  of  propagation  being  less  than  the  free 
space  velocity  and  from  an  increase  in  path  length  brought  about  by  the  refractive  bending  of  the  ray, 
reduces  to 


AR=  1  (Rj*nj  )  -  Rq^ 

7=0 


(38) 


where  the  distance  Rj  is  given  by 


^7'  =  +  'f+l  ~  cos  Qj .  (39) 

Since  the  validity  of  Snell’s  law  or  Bouger’s  law  depends  upon  the  thickness  of  layers,  the  performance 
of  the  stratified  model  largely  depends  on  the  number  of  layers  in  the  process.  The  total  number  of  layers 
is  45  in  this  report  and  is  specified  in  the  following  manner; 

/i=0-100m  10  layers  with  10-m  interval 

/i  =  100  -  1,000  m  9  layers  with  100-m  interval 

h  =  1,000  m  and  above  26  layers  with  1,000-m  interval.  (40) 

The  main  reason  for  this  division  is  based  on  the  fact  that  the  refractive  effect  is  small  above  10  km  from 
the  surface  of  the  Earth.  The  test  has  been  performed  in  order  to  validate  this  argument  of  45  layers  by 
computing  range  and  angle  errors  on  different  environments  with  spatial,  temporal,  and  geographical 
variations.  This  approach  provides  both  range  and  angle  errors  with  straightforward  mathematical  and 
computational  expressions. 
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3.2.3.  Goad  Model  [6] 

This  model  results  from  the  combination  of  the  Hopfield  [2]  and  Saastamoinen  [19]  models.  Goad 
modified  the  value  of  the  tropospheric  height  into  a  Taylor’s  series  approximation  in  terms  of  range  rather 
than  a  quartic  form  as  proposed  by  Hopfield  as 

/i  =  r  sin  0  +  cos^  >  (41) 

where  h  is  height,  r  the  range,  0  the  elevation  angle,  and  Og  the  semimajor  axis  of  the  Earth. 

The  range  correction  AR  is  computed  as: 

AR  =  10-6  Jq No (r)  +  10-6  Ni{r)dr 

=  10-^  Ni  (0)  [flij  ri  +  (a2j  /  2)  r^  + ...  +  (ag j  /  9)  ,  (42) 

where  Nq  is  the  surface  dry  refractivity,  and  the  wet  refractivity 

r,  =  -^(og cos^  sinO,  for  7  =  0,1 

hQ  =  5(0.002277)/>/A/o*10-6 

Til  =  [5(0.077) /(iVi  *10~^)]*  (1255/7’+ 0.5}*e 

ai/=l,  0.21=4  0, 

o^I=6ai  2  +4bi  a^j  =  4  a,-  (a?  +  3bi ) 

O51  =  af  +12af  bi  +6bf 

«6/  =4a,-  bi  (af  +3/7,  ) ,  =  bf  {6af  +4/?;) 

Ogi  =40,  bf ,  Ogi  =  bf 


Oi  =  sin  Q/hi, 


e  =  6.108*7577  *exp[(17. 15  7- 4684)7(7  -  38.45)1 
7  =  0:  dry  refractive  component,  7=1:  wet  refractive  component 


14 


Junho  Choi 


Ar:  range  correction  in  m,  subtract  from  pseudo  range  or  carrier  phase 

T:  surface  temperature  in  K  (  =  °C  +  273.16) 

p:  atmospheric  pressure  in  mbar 

e:  water  vapor  partial  pressure  in  mbar 

0:  elevation  angle  tangent  to  the  horizon 

a^:  semimajor  axis  of  the  Earth  ellipsoid 

FH :  relative  humidity  as  a  fraction  of  1 .0 
Similarly  the  elevation  angle  correction  is  obtained  as 

A0(r)  =  jo^i0  =  4coseo  I  {XuR  +  X2i  12  + ...  + X-j  i  /l]  ,  (43) 

i=0  J 


where 

^u  =  i 
■^2,1  = 

•^3,1  =  3(a,-2  +  bf) 

X^j  =  flj'  (6  bi  +  ap-) 

■^5,1=3  bi^bi  +  aP') 

X6,,  =  3a,fc,2 

Xij  =  bi^ 

and  other  parameters  are  already  defined  above.  As  pointed  out  before,  this  model  covers  both  below  5° 
elevation  angle  and  between  5°  and  20°  in  contrast  with  Hopfield’s  model,  which  covers  only  above  5° 
elevation  angles.  Test  and  evaluation  of  this  model  was  somewhat  limited  to  a  couple  of  data  samples. 

3.2.4.  Blake  Model  [5] 

Blake  developed  a  technique  to  find  the  position  coordinates  of  the  point  on  the  Earth’s  surface  that 
lies  directly  below  the  target,  which  amounts  to  finding  the  ground  range  (distance  from  the  antenna  to 
the  target  point  measured  along  the  Earth’s  surface  at  sea  level).  Further  quantities  of  interest  are  the 
straight  line  distance  from  antenna  to  target  and  the  true  elevation  angle  (as  opposed  to  the  apparent  angle 
indicated  by  the  antenna).  Consider  the  ray  tracing  equation,  which  is  Snell’s  law  for  a  spherically 
symmetric  medium. 


cos  0  =  [«o  cos  00  ]/n(/r)  (l  +  /  tq  )  ,  (44) 

where  0o  is  the  initial  elevation  angle  at  /i  =  0,  tq  the  distance  from  the  Earth’s  center  to  the  initial  point, 
no  the  refractive  index  at  /i  =  0  height,  and  h  the  ray  height  above  its  initial  point.  Substituting  Eq.  (44) 
into  Eq.  (24),  one  can  obtain  the  range  error  as 
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R  (^i,6o)  =  j  l--^|no  coseo/[«(l  +  /i/ro)]} 


(45) 


A  similar  procedure  yields  for  the  ground  range  as 

G  {hi,Qo)=  +  +  COS0O  ]  ^  -1  j  .  (46) 

In  Fig.  2,  the  quantity  R'  is  the  geometric  length  of  the  ray  path.  The  corresponding  distance  R,  measured 
by  a  radar  that  assumes  the  wave-propagation  speed  to  be  that  for  free  space  (the  speed  of  light 
c  =  2.997925  x  10*  m/s),  is  related  to  R'  by  the  differential  expression 

dR  =  ndR'.  (47) 


Fig.  2  —  Ray  path  geometry  in  an  atmosphere  spherically  symmetric  with  respect  to  Earth’s  center 


The  quantity  R  given  by  Eq.  (45)  is  the  radio  path  range.  The  distance  G  given  by  Eq.  (46)  is  the  ground 
range  equal  to  The  particular  refractive  index  profile  for  which  computations  have  been  carried  out  is 

of  the  form 


n(h')  =  1  +  exp  -  /ij)] , 


(48) 
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where  kisa  decay  constant,  and  (l  +  p^)  is  the  value  of  n  at  h'  =  h^\  that  is,  is  a  reference  height.  If 
the  ray  origin  is  at  some  height  (antenna  height)  and  h  is  measured  with  respect  to  the  antenna,  the 
expression  (Eq.  (48))  becomes 


n(h)  =  1  +  po  e\p(-kh) ,  (49) 

where  Po  =  p^  exp  [-k{h^  -  /i,)]  .  (50) 

At  small  values  of  Bq,  the  integrands  of  Eqs.  (45)  and  (46)  are  subject  to  the  loss  of  accuracy  in  the 
vicinity  of  h  =  0.  Blake  manipulated  the  integral  equations  for  numerical  computations  as 

R  (/ii,6o)  =  Ja=o  (1+  h/rQ)dhy^fu  +  V  +  W  +  VW,  (51) 

G(/ii,0o)  =  J^=o  {(1+  Po)(cos0o)fi?^}/[(1+  /i/ro)(l+  /i/ro)VM  + v  +  w  +  vw,]  (52) 

where  n  =  1  +  Pq  exp(-A/j) 

«=  (l  +  Po)^sin2  0o-2po-  Po 
V  =  Po  exp  i-kh)  +  Po  e\p(-2kh) 
v/  =  2h/rQ  +  h^/ri. 


These  formulations  are  obtained  by  replacing  cos^  0q  with  its  equivalent,  ^1  -  sin^  0o  j.  At  0o  =  0,  the 
quantity  «  is  a  small  negative  number,  but  above  about  0o  =  1.5°,  it  becomes  positive.  In  the  region  close 
to  00  =  0  and  h  =  0,a  special  procedure  is  required.  As  pointed  out,  Blake’s  model  took  into  consideration 

the  low  elevation  angle  (below  5°)  to  compensate  the  shortfalls  of  Hopfield’s  model,  which  is  valid  only 
above  5  °.  Blake  also  developed  a  technique  to  compensate  the  range  error  for  the  over-the-horizon 
atmospheric  effects  by  using  Eq.  (46)  rather  than  Eq.  (45)  alone  (note:  most  models  use  a  single 
expression  for  the  range  integration). 

Similarly,  equations  for  the  true  geometric  range  /?,,  and  the  true  elevation  angle  0,  are  manipulated 
as 


I?,  =  V  ^hi  +4ro  (ro  -t-/zi)sin^  (G/2/'o)  , 
0j  =  sin  I  Rf  +  /  2rQ  R^  —  R^  /  2/^  j 

0,  =  7t/2  -  sin“^  ||(ro  +  /ii )  sin  (G  /  tq  )}  / 


(53) 

for  00  <  n/4, 

(54) 

for  00  >  71/4. 

(55) 

Note  that  the  accuracy  of  all  of  these  results  is  heavily  dependent  on  the  accuracy  of  the  calculation  of  G. 
The  elevation  angle  error  can  be  obtained  readily  by  subtracting  the  true  elevation  angle  of  Eqs.  (54)  and 
(55)  from  the  actual  elevation  angle. 
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3.2.5.  Cain’s  Model  [20] 


This  model  was  developed  by  D.L.  Cain  of  the  Jet  Propulsion  Laboratory  (JPL),  Pasadena, 
California,  in  the  late  1960s  or  early  1970s,  for  deep-space  communication-link  correction.  A  direct 
reference  for  Cain’s  work  could  not  be  found,  but  a  brief  outline  of  his  model  was  given  in  Moyer  [20] 
and  Gallini  [21],  This  model  is  one  of  the  closest  models  to  the  modified  exponential  model.  Thus,  this 
model  is  included  with  other  model  performances  for  comparison  purposes.  The  model  has  the  form  for 
range  delay  AR  as: 


AR=  (C4/340*Af^)/[(C2sin£-i-Ci)*sin£+Cop  [m] 


(56) 


=  station  surface  refractivity  index 


E  =  geometric  elevation  angle 
Cq  =  0.06483  [unitless] 

Cj  =1.0  [unitless] 

C2  =  0.0  [unitless] 

C3  =  1 .4  [unitless] 

C4  =  0.0018958  [km]  . 


This  model  has  been  implemented  in  the  Air  Force  Satellite  Control  Network  Tracking  and  Orbit 
Determination  program  and  uses  a  table  of  monthly  averages  for  the  surface  refractivities.  According  to 
Ref.  21,  current  monthly  averages  are  somewhat  low,  and  the  formula  itself  results  in  inaccurate 
refractivities.  For  a  given  refractivity,  the  formula  produces  a  correction  of  14  m  at  10°  elevation  angle, 
which  is  a  reasonable  value. 


3.2.6.  Case  1  Model 

Since  no  reference  for  this  model  could  be  found,  we  called  it  Case  1  Model.  This  model  has  been 
adopted  in  the  program  for  more  than  a  decade  without  any  proof  or  verification  of  its  validity.  The  origin 
of  this  model  was  derived  from  Figs.  1-4  and  1-5  of  Millman  [4]  for  both  elevation  angle  and  range  error 
plots,  respectively,  by  segmenting  into  three  parts  and  then  computing  coefficients  through  least-square- 
fit  approach.  This  is  similar  to  Cain’s  model  except  for  separation  into  three  parts  based  on  elevation 
angle.  The  range  delay  AR  in  meters  is  given  by 


AR  =  0.3048006/10.003589  +  (0.087605  sin  E  +  0.19696793  sin  2  E)]  for  E  <  5°, 

AR  =  0.3048006/[0.002129  +  0.12158  sin  R]  for  5  <  E  <  30°, 

AR  =  0.3048006/[0.03  +  0.08  sin  E]  for  R  >  30° , 


(57a) 

(57b) 

(57c) 
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where  E  denotes  the  elevation  angle,  and  c  is  the  speed  of  light,  299.792458  *10^  km/s.  The  time  delay 
can  be  calculated  readily  by  dividing  range  error  by  c.  Note  here  that  both  this  model  and  Cain’s  model 
do  not  provide  an  angle  error  formula  and  are  similar. 

3.3.  Proposed  Approach 

Here  we  propose  the  approach  that  combines  both  the  modified  exponential  refractivity  model  and 
the  stratified  range/angle  error  model  presented  above.  This  approach  gives  more  accurate  range  and 
angle  error  performance  in  the  low  elevation  angle  below  10°  in  comparison  with  that  of  other  models. 
Refractivity  can  be  generated  without  any  limitation  of  altitude  from  ground  to  the  top  of  the  tropospheric 
layer  and  is  closer  to  the  empirical  data  than  any  other  approach.  This  model  may  also  be  implemented 
easily  for  a  combination  of  tropospheric  and  ionospheric  effects  and  can  be  extended  to  the  ionospheric 
layer  without  difficulty.  The  total  storage  (only  7  Mbps)  and  programming  requirements  (<  100  lines  of 
coding)  are  minimal  except  for  the  database  of  the  2.5°  x  2.5°  worldwide  reference  heights,  which  can  be 
updated  every  5  to  10  yr  depending  upon  the  user’s  accuracy  requirements.  Test  and  analysis  results  are 
given  in  the  next  section. 

4.  PERFORMANCE  COMPARISON  OF  EACH  MODEL 

Evaluation  of  model  performance  is  emphasized  mainly  on  temperature,  relative  humidity, 
refractivity  distribution  and  its  gradient,  time  delay,  range  error,  and  elevation  angle  error  over  more  than 
130  areas  of  interest  (AOIs)  with  worldwide  coverage  in  geographical,  climatological,  surface  altitude, 
seasonal,  diurnal,  hourly,  marine,  and  polar  regions.  Methodology  of  data  representation  varies  with 
characteristics  of  data  contents  by  graphs,  tables,  contours,  and  colors.  Since  the  volume  of  analysis 
results  are  enormous,  the  majority  of  data  is  included  in  the  Appendixes.  Data  are  divided  into  three  parts: 
10  AOIs,  46  AOIs,  and  130  AOIs,  depending  upon  the  number  of  variables  for  the  comparison  purposes. 
The  comparison  is  divided  into  two  parts — the  first  on  refractivity  and  meteorological  parameters 
(temperature,  relative  humidity,  pressure,  and  refractivity  and  the  second  on  propagation  errors.  Figure  3 
shows  the  10  AOI  locations  selected  for  simplification  of  data  handling  to  analyze  first-hand  geographical 
performance  variations.  Figure  4  presents  46  AOIs  selected  to  study  further  details  and  trends  of  model 
performance  for  different  geographical  and  climatological  regions.  Table  1  presents  the  total  number  of 
AOIs  examined  in  this  report  for  evaluation  of  model  performance  and  data  quality  verification  and 
validation.  Table  1  contains  full  descriptions  and  abbreviations  of  all  AOIs  considered  here  with  latitude 
and  longitude  boundaries  and  number  of  grids  in  that  area.  Note  here  that  some  portions  of  data  have  been 
interpolated  or  extrapolated  based  on  the  availability  of  raw  data  because  original  data  have  portions  of 
missing  or  erroneous  data  for  some  grids  or  areas.  All  data  analyses  presented  in  this  report  are  based  on 
three  categories  of  data  groups:  10,  46,  or  130  AOIs. 

4,1  Refractivity  and  Meteorological  Parameters 

Figures  5  and  6  show  the  worldwide  ECM  data  distribution  of  refractivities  for  Febraary  and  August, 
respectively,  for  15°  intervals  of  latitude  and  45°  intervals  of  longitude.  Higher  refractivity  regions  for 
August  (Fig.  5)  are  in  the  Northern  Hemisphere  with  higher  refractivity  regions  in  the  Southern 
Hemisphere  during  Febraary  (Fig.  6).  Refractivities  for  areas  of  30°  latitude  of  Northern  Hemisphere  and 
30°  latitude  of  Southern  Hemisphere  are  much  higher  than  the  worldwide  average,  332  to  340  N-units 
depending  on  time  of  the  day.  Table  2  presents  46  AOIs  with  latitude,  longitude,  mountain  (M),  land  (L), 
water  (W),  continents,  and  climates.  These  46  AOIs  are  selected  to  cover  weighted  areas  of  desert,  ocean, 
rain  forest,  tropical,  temperate,  and  polar  regions  for  fair  evaluation  of  model  performance  equally  and 
evenly.  Table  3  shows  the  monthly  surface  refractivity  for  46  AOIs  of  ECM  surface  data.  Note  here  that 
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Table  1  —  Selected  Areas  of  Interest 
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Table  1  (Continued)  —  Selected  Areas  of  Interest 


67 

Greater  Australia  Continent 

GRAUS 

-30.0 

-10.0 

110.0 

160.0 

160 

68 

Northern  Australia:  Tanami  Desert 

NAUS 

-25.0 

-15.0 

130.0 

137.5 

12 

69 

Antarctic  Circle 

ANT 

-90.0 

-65.0 

0.0 

360.0 

1440 

70 

Northern  China  Desert 

NCHD 

35.0 

47.5 

75.0 

95.0 

40 

71 

West  Indian  Coast 

WIC 

■■1X1 

25.0 

60.0 

75.0 

36 

Latitude  (deg) 
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Refractivity  (N) 

Fig.  6  —  ECM  worldwide  refractivities  (February) 


Table  2  —  Definitions  and  Characteristics  for  46  Areas  of  Interest 
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Table  3  —  Monthly  Surface  Refractivity  (N)  for  46  Areas  of  Intrest  from  ECM  Surface  Data 
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the  surface  refractivity  is  derived  from  extrapolation  of  the  original  raw  data  supplied  by  the  originator. 
Table  4  shows  monthly  temperature  distribution  for  46  AOIs  of  ECM  surface  data.  Tables  5  and  6  present 
monthly  surface  relative  humidity  and  surface  pressure,  respectively,  for  46  AOIs  of  ECM  surface  data. 

Figures  7  and  8  show  the  ECM  surface  relative  humidity  worldwide  contour  map  monthly  10-year 
averages  for  February  and  August,  respectively,  to  present  the  contrast  between  the  coldest  and  hottest 
months  of  the  year.  Figures  9  and  10  present  ECM  surface  refractivity  worldwide  contour  maps  for 
February  and  August  monthly  10-year  averages,  respectively.  High  refractivity  areas  slowly  move  to  the 
Northern  Hemisphere  with  low  refractivity  areas  noticeably  in  the  Northern  Hemisphere.  Figures  1 1  and 
12  present  HIRAS  surface  average  refractivity  worldwide  contour  maps  for  February  and  August, 
respectively.  These  HIRAS  contour  maps  noticeably  contrast  with  those  of  Figs.  9  and  10.  ECM  data  are 
generated  from  10-year  monthly  average  climatology  data  while  HIRAS  data  are  generated  from  six 
hourly  9-year  average  empirical  data. 

Figure  13  through  18  present  direct  interrelationship  among  temperature,  relative  humidity,  and 
refractivity  for  the  months  of  February  and  August  with  three  AOIs.  Figures  13  and  14  show  the 
interrelationship  among  meteorological  parameters  in  the  Washington,  D.C.  area  for  the  month  of 
February  to  represents  winter  season  and  for  the  month  of  August  to  represent  summer  season.  Diurnal 
variations  of  both  temperature  and  relative  humidity  are  more  dynamic  in  February  than  in  August.  This 
implies  more  variations  of  refractivity  in  February  than  in  August.  Figures  15  and  16  show  the 
interrelationship  among  meteorological  parameters  in  the  neighborhood  of  the  Aleutian  Islands  of  Alaska 
(NAK).  Diurnal  variations  of  temperature  and  relative  humidity  are  similar  to  those  in  the  Washington, 
D.C.  area  for  both  February  and  August.  Figures  17  and  18  show  the  interrelationship  among 
meteorological  parameters  for  the  Amazon  rain  forest  (AMFOR)  area.  Diurnal  variations  of  both 
temperature  and  relative  humidity  are  similar  for  the  months  of  both  February  and  August.  Details  for 
other  areas  are  included  in  Appendixes  A  to  D. 

4.2.  Range  or  Time  Delay  and  Angle-of-arrival  Errors 

Time  delays  and  angle-of-arrival  errors  based  on  empirical  databases  and  selected  models  appear  in 
Tables  7  through  11.  Table  7  provides  time  delays  for  HIRAS  empirical  data  and  six  other  models  from 
the  horizon  to  the  elevation  angle  of  20°.  The  variation  of  time  delays  on  seven  different  cases  is 
noticeable.  Results  of  Table  7  indicate  that  the  choice  of  model  dominates  the  performance  of  range 
errors  and  time  delays.  Table  7  also  shows  that  time  delay  differences  are  minimal  above  the  20° 
elevation  angle  in  comparison  with  those  of  lower  elevation  angles.  Table  12  tabulates  characteristics  of 
selected  models  with  given  weather  and  geographical  parameters  and  newly  calculated  refractivity  related 
parameters.  Some  models  do  not  provide  all  three  of  the  outputs  of  time  delay,  range  error,  and  angle-of- 
arrival  error.  Figures  19  through  22  present  time  delays  and  angle-of-arrival  errors  for  different  models 
with  empirical  data  of  February,  May,  August,  and  November  in  the  Amazon  rain  forest.  Time  delay 
performance  of  the  exponential  model  is  close  to  the  empirical  data  (stratified  model)  within  2%  of  sum- 
squared  errors  with  empirical  data  while  Hopfield,  Goad,  and  Cain’s  model  performances  are  worse  with 
20%  to  30%  sum  squared  errors;  Blake  model’s  performance  is  within  10%  sum  squared  error.  Angle-of- 
arrival  error  for  the  exponential  model  is  worse  than  that  of  stratified  model  below  1°  elevation  angle  with 
20%  while  the  angle  error  of  the  Goad  and  Blake  models  is  worse  than  that  of  the  stratified  model  with 
70%  sum  squared  error  below  1°  elevation  angle.  It  is  noted  in  these  figures  that  both  time  delays  and 
angle-of-arrival  errors  are  at  least  50  ns  higher  than  other  areas  all  year  around  for  time  delays  and  0.2°- 
0.3°  elevation  angle  errors  when  compared  with  other  worldwide  areas.  These  imply  that  RF  wave 
bending  or  time  delay  depends  highly  on  both  temperature  and  humidity  since  both  temperature  and 
humidity  are  always  high  in  the  rain  forest.  Figures  for  other  AOIs  are  included  in  Appendixes  E  to  J  for  a 
complete  picture  of  the  different  geolocation  and  climatology. 
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Performance  Comparison  of  Tropospheric  Propagation  Models 


New  Caiedonia  (NEWCAL)  (43) _ 1008.93  1007.90  1010.31  1012.37 

New  Zealand  (NEWZEA)  (44) _ 1010.14  1012.65  1014.96  1017.07 

Rio  de  Janerlo  (RIODEJ)  (45) _ 1010.69  1011.67  1012.21  1014.51 

Falklands  (FALKLD)  {46)  999.0ol  1000.00)  1002.051  1001.^ 
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Refractivity  Temperature  Humidity  Refractivity  Temperature  Humidity 
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Time  Ranges  (Days) 
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Fig.  15  —  MFR  time  series  surface  data  -  NAK  -  August  1995 
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Pig,  17  —  MFR  time  series  surface  data  -  AMFOR  -  February  1995 
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Fig.  18  —  MFR  time  series  surface  data  -  AMFOR  -  August  1995 
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Table  7  —  Time  Delay  for  D.C.  Area  using  HIRAS  Data  and  Various  Models 

TIME  DELAY:  HIRAS  Data  DC  Araa  Juna  OOOOKra 
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Table  8  —  Time  Delay  (ns)  for  Selected  Areas  of  Interest 
MRF,  Goad,  and  Exponential  Model  for  15  February  1995 
(0000,  0600,  1200,  and  1800  h) 
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Table  9  —  Time  Delay  (ns)  for  Selected  Areas  of  Interest 
MRF,  Goad,  and  Exponential  Model  for  15  August  1995 
(0000,  0600, 1200,  and  1800  h) 
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Table  10  —  Angle  Error  (degrees)  for  Selected  Areas  of  Interest 
MRF,  Goad,  and  Exponential  Model  for  15  February  1997 
(0000,  0600, 1200,  and  1800  h) 
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Table  11  —  Angle  Error  (degrees)  for  Selected  Areas  of  Interest 
MRF,  Goad,  and  Exponential  Model  for  15  August  1997 
(0000,  0600,  1200,  and  1800  h) 
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Table  12  —  Parametric  Requirements  for  Tropospheric  Models 
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X 
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Elevation  Angle 

X 

X 

X 

X 

X 

X 

Pressure  (mb) 
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X 
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X 
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|Wet  Height  of  Tropospheric 

X 

Calc’d 
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Fig.  19  —  ECM  database  -  AMFOR  -  February 


Angle  Error  (Degrees)  Time  Delay  (Nanoseconds) 
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Angle  Error  (Degrees)  Time  Delay  (Nanoseconds) 
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4,3.  Azimuth  Angle  Dependence  of  RF  Ray  Bending 

It  has  been  assumed  that  the  refractive  index  of  the  atmosphere  in  ray-tracing  studies  is  spherically 
stratified  with  respect  to  the  surface  of  the  Earth.  Thus,  the  effect  of  refractive  index  changes  in  the 
horizontal  direction  is  normally  not  considered.  Equations  (27)  through  (30)  are  subject  to  the  following 
assumptions  of  ray  tracing: 

(1)  The  refractive  index  should  not  change  appreciably  in  a  wavelength. 

(2)  The  fractional  change  in  the  spacing  between  neighboring  rays  (initially  parallel)  must  be  small 
in  a  wavelength. 

(3)  The  refractive  index  stracture  is  horizontally  homogeneous. 

(4)  The  refractive  index  is  a  function  only  of  height  above  the  surface  of  a  smooth  and  spherical 
Earth. 

Neglecting  the  effect  of  horizontal  gradients  seems  reasonable  in  the  tropospheric  region  because  of 
the  relatively  slow  horizontal  change  of  refractive  index  in  contrast  to  the  rapid  decrease  with  height. 
Bean  and  Dutton  [1]  performed  two  experiments  (one  over  the  Canterbury  in  New  Zealand  and  the  other 
at  Cape  Canaveral,  Florida)  to  investigate  assumption  (3)  (above)  for  horizontal  changes  of  the  refractive 
index.  Their  conclusion  was  not  clear  and  mixed  with  the  emphasis  of  ducting  studies.  Vogel  [22] 
investigated  monthly  variations  of  refractivities  and  ducting  for  horizontal  effects  of  radio  wave 
propagation. 

The  importance  here  of  this  azimuth-angle  dependence  is  to  achieve  a  high  accuracy  in  the  low 
elevation  angle  (or  horizon)  for  some  over-the-horizon  radio  communication  applications.  For  example, 
the  search  and  surveillance  mission  over  the  horizon  always  encounters  the  atmospheric  effect 
(tropospheric  region)  on  RF  propagation  over  an  entire  azimuth  direction  (over  360°).  If  one  applies  one 
uniform  refractive  index  over  entire  azimuth  direction,  one  can  obtain  erroneous  results  over  each 
different  azimuth-angle  measurement.  Figures  23  to  26  show  time  delay  plots  over  each  90°  azimuth 
direction  over  360°  from  the  months  of  February  to  November  for  six  hourly  averages  between  1988  and 
1994  in  the  Washington,  D.C.  area.  The  patterns  of  time  delays  over  the  entire  azimuth  direction  vary 
hourly  and  seasonally.  The  variation  of  time  delays  for  the  months  of  May  and  August  is  sharply 
distinguishable  from  other  monthly  and  hourly  observations.  The  analysis  is  performed  over  several 
regions  to  find  conclusive  evidence.  If  the  elevation  angle  is  above  5°,  the  azimuth-angle  dependence  is 
almost  negligible.  If  the  elevation  angle  is  between  0°  and  3°,  both  time  delay  and  range  error  should  be 
carefully  calibrated  or  calculated  for  correct  tracking  and  navigation  applications,  including  geolocation 
accuracy.  Additional  graphs  are  included  in  Appendix  K  for  Teheran,  Iran,  and  Ahaggar,  Algeria. 

5.  CONCLUSIONS  AND  RECOMMENDATIONS 

The  major  factors  of  model  comparison  criteria  are  based  on  the  accuracy,  minimum  level  of 
database  or  information  requirement,  and  feasibility  of  real-time  data  application.  Each  model  has 
advantages  and  disadvantages  in  each  category.  As  we  pointed  out  in  Section  4,  most  of  the  comparisons 
are  performed  through  time  delay,  range  error,  and  angle-of-arrival  error  for  different  AOIs  on  different 


TIitw  Delay  (ns)  Time  Delay  (ns) 


Performance  Comparison  of  Tropospheric  Propagation  Models 


Time  Delay  (ns) 


0.0  50.0  100.0  150.0  200.0  250.0  300.0 

Azimuth  (deg) 

Fig.  25  —  D.C.,  August  (HIRAS  surface  data — 0°  elevation  angle) 


Fig.  26  —  D.C.,  November  (HIRAS  surface  data — 0°  elevation  angle) 
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climatology  to  induce  reasonable  conclusions.  A  modified  exponential  model  is  proposed  here  as  the  best 
performer  among  many  models  evaluated  based  upon  the  three  criteria  level  of  accuracy,  minimum  level 
of  database  including  surface  weather  data,  and  real-time  data  applicability.  The  modified  exponential 
model  outperforms  in  the  accuracy  improvement  of  range  and  time  delay  errors  while  surface  weather 
data  and  reference  height  are  required  only  to  run  ray-tracing  algorithms  for  range,  time  delay,  and  angle- 
of-arrival  error.  This  model  also  instantaneously  accepts  real-time  weather  data  anytime  and  anywhere  in 
the  world.  Several  technical  evaluation  efforts  are  currently  under  way  in  programs  and  projects  in  many 
agencies  and  will  be  extended  further.  Preliminary  results  show  promise  for  implementation  to  real 
operating  systems  in  both  tracking,  surveillance,  and  navigation  applications. 
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Appendix  A 

REFERENCE  (SCALE)  HEIGHTS  BY  LATITUDE, 
LONGITUDE,  AND  MONTHS 


Latitude  is  semented  by  15°  from  -60°  to  +60°and  by  30°  from  60°  to  90°.  Longitude  is  divided  by 
90°  from  -180°  to  +180°  for  entire  globe. 
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Reference  Height  (m) 


54 


Junho  Choi 


I  ?  Ill 

o>  o  5, 
o  a> 

O  Q  & 

m  at* 


?  5 
O  LU  O 
o»  o  “ 


S  °  b| 


^  >  LU 

.  fmh 

a  o  ^  CD 

o  o  O 


^5  ' 

_  r  uj  . 

*  o  ^  . 
o  o> 

ffz"' 

o  o  6 


*o 

3 

I 


S 


is 


eo 

■M- 

y- 

cc 

h- 

CO 

CM 

Oi 

o 

CM 

CD 

CD 

r*. 

o 

CD 

o 

CM 

CO 

in 

CO 

cd 

cr 

cr 

(X 

c 

IT 

cr 

CV 

IT 

fv 

a 

tf 

IT 

CC 

cc 

in 

CO 

CO 

CO 

CO 

CD 

CO 

CO 

CO 

f- 

Oi 

CO 

o 

a> 

o 

in 

CD 

o 

CO 

O) 

o 

CM 

CM 

CM 

CO 

CM 

in 

in 

CM 

ir 

CC 

■««“ 

CV 

in 

CO 

IT 

C£ 

IX 

M 

lO 

cc 

in 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

o 

o 

CO 

in 

CO 

T* 

o 

o 

CO 

CD 

CO 

CD 

CO 

CO 

eo 

CM 

m 

CO 

N. 

03 

Oi 

c 

o 

CV 

in 

CV 

O) 

IT) 

in 

in 

CO 

CO 

CO 

CO 

CO 

CO 

CD 

CO 

o 

■M- 

O) 

CD 

■M- 

M- 

o> 

'«■ 

CO 

'M' 

CO 

CO 

in 

in 

in 

to 

CO 

CO 

o 

in 

o 

a 

to 

m 

m 

CO 

in 

CD 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

T- 

CO 

y- 

O) 

in 

fv 

CO 

in 

in 

CO 

CO 

CM 

CO 

in 

o 

o 

CO 

o 

CM 

o 

CD 

K 

CO 

CM 

CO 

a> 

CM 

CO 

in 

in 

CO 

in 

in 

K 

a> 

CO 

o> 

CO 

CO 

CO 

CO 

CO 

CO 

T- 

00 

in 

CM 

o 

CO 

▼- 

CD 

CM 

CO 

CO 

cd 

to 

■M- 

o 

CO 

d 

CO 

CO 

CO 

in 

CD 

in 

in 

in 

CO 

in 

in 

f»- 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CD 

CO 

CM 

CO 

CO 

CO 

to 

in 

CM 

-M- 

CM 

to 

o 

m 

h* 

in 

in 

CM 

CT 

m 

T-p 

CM 

CO 

o 

Oi 

CD 

CM 

CO 

CO 

CO 

CO 

O) 

■V 

in 

CD 

CO 

to 

in 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

CM 

o> 

t- 

o 

CM 

o 

CM 

CM 

CO 

r“ 

r«» 

CM 

to 

in 

CD 

o 

CM 

O 

o 

CD 

in 

■M- 

h- 

CM 

-M- 

CM 

CM 

CD 

’♦ 

CD 

CO 

'*■ 

'M- 

in 

in 

CD 

CO 

CO 

CO 

CO 

CO 

CO 

CO 

*« 

m 

"ct 

“ 

aJ 

To 

:  CD 

in 

O) 

CO 

CM 

CO 

o 

CO 

CD 

ed 

CD 

in 

CM 

d 

cd 

in 

in 

'  in 

Q 

o 

in 

1  eo 

CO 

r«. 

1  CO 

'M- 

in 

CD 

CO 

CO 

'  CO 

CO 

CO 

to 

to 

To 

To 

■M' 

To 

— 

"d 

h- 

CO 

CD 

CD 

O) 

CM 

CD 

cd 

in 

cd 

d 

V 

V 

cd 

CO 

CO 

eo 

ID 

'  CM 

o 

m 

CM 

in 

r««. 

o 

■M- 

in 

CO 

in 

in 

CO 

CD 

r«* 

CO 

CO 

to 

CD 

to 

To 

"o 

To 

To 

Tm 

"S 

in 

CO 

CM 

CO 

CO 

CM 

o 

CD 

d 

CD 

K 

d 

in 

eo 

in 

in 

CM 

o 

CO 

in 

to 

m 

m 

CO 

CO 

in 

CD 

CO 

“ 

h- 

CO 

to 

CD 

CD 

to 

“cD 

To 

CD 

Oi 

CO 

o 

CC) 

o 

eo 

in 

CM 

K 

CM 

d 

CD* 

in 

eo 

CD 

CO 

in 

CD 

in 

in 

CO 

CO 

in 

CO 

CO 

CD 

CD 

CO 

CO 

CO 

f 

o 

? 

lii 

UJ 

h 

f 

d 

f 

UJ 

UJ 

d 

CD 

O 

h 

a) 

CD 

o 

d 

to 

O 

o 

CD 

O 

o 

CD 

? 

$ 

O 

o 

UJ 

5 

i 

o 

o 

UJ 

o 

€0 

h 

CD 

d 

d 

CD 

d 

ed 

d 

CD 

d 

d 

a 

z 

d 

& 

o 

o 

« 

o 

< 

S 

I 

s 


.r> 

■M- 

03 

CD 

Jx 

a 

a 

a 

a 

a 

a 

a 

a 

a 

CM 

(D 

CO 

r«* 

CO 

CM 

CM 

o 

CM 

o 

CM 

a 

cc 

o 

o 

CM 

CO 

CO 

■M- 

o 

Oi 

tx 

V- 

a 

CM 

CM 

o 

a 

■vr 

CM 

a 

CD 

a 

cr 

u 

x) 

rx 

<* 

CV 

cr 

u 

cc 

C 

O 

a 

QC 

c 

■M 

in 

in 

o 

C 

c 

o 

c 

to 

CO 

CO 

fx- 

rx. 

rx 

|X 

|X 

rx 

a 

a 

IX 

a 

a 

a 

a 

CO 

CD 

Oi 

CO 

T- 

CM 

a 

a 

rx 

a 

y- 

a 

a 

rx 

a 

« 

in 

CM 

CD 

<D 

O 

a 

CM 

a 

a 

a 

o 

a 

CD 

CD 

o 

CO 

a 

'M' 

o 

CM 

rx 

CM 

CM 

CO 

o 

cd 

(X 

CX 

m 

rx 

o 

in 

•M" 

|X 

o 

o 

c 

a 

CX 

a. 

CX 

CV 

in 

fC 

in 

CC 

O! 

C 

C 

o 

c 

o 

CO 

CO 

CO 

rx. 

rx 

rx 

rx 

rx 

rx 

a 

a 

rx 

a 

a 

a 

a 

CM 

in 

OI 

CO 

CO 

CO 

"If 

o 

a 

a 

|X 

a 

CM 

u> 

o 

OI 

CM 

a 

a 

CO 

o 

fx 

a 

a 

a 

rx 

'W- 

■M- 

•r- 

a 

o 

rx 

a 

CM 

a 

CM 

rx 

CM 

a 

a 

T“ 

CM 

in 

cc 

CM 

■M- 

cc 

CV 

a 

cr 

a 

ir 

in 

a 

a 

cr 

cc 

rx 

a 

rx 

c 

c 

c 

c 

c 

o 

a 

CO 

K 

rv 

rx. 

rx 

rx 

rx 

rx 

a 

a 

a 

a 

a 

a 

rx 

rx 

CM 

CO 

to 

CO 

CO 

o 

Oi 

a 

a 

a 

a 

CM 

rx 

rx 

a 

•M- 

CD 

CD 

in 

CD 

Ol 

a 

a 

a 

■M' 

a 

•*- 

•r* 

a 

a 

a 

O 

CM 

T- 

CO 

CO 

CO 

a 

a 

a 

a 

a 

a 

fx 

a 

cd 

■M” 

03 

CV 

CM 

a 

a 

a 

a 

CT 

CM 

a 

a 

a 

a 

fx 

O 

CM 

rx 

rx 

io 

fx 

a 

c 

a 

a 

a 

a 

r. 

r*. 

rx 

rx 

rx 

rx 

a 

to 

a 

a 

rx 

rx 

IX 

rx 

T- 

CO 

'M' 

fx 

rx 

rx 

EO 

a 

a 

CM 

o 

CM 

O 

a 

Oi 

in 

in 

CO 

CM 

CO 

CM 

o 

CT 

a 

a 

a 

a 

a 

in 

CM 

rx. 

CO 

CO 

0) 

a 

a 

CD 

"«• 

■Vf 

a 

d 

CM 

CD 

o 

CM 

a 

O 

a 

a 

CM 

o 

CO 

a 

rx 

fx 

O 

CM 

■M- 

rx 

rx 

rx 

a 

o 

a 

a 

a 

a 

r». 

rx. 

rx 

rx 

rx 

rx 

to 

a 

a 

a 

rx 

rx 

rx 

fx 

CO 

CD 

CM 

CO 

a 

T- 

a 

a 

n 

rx 

a 

a 

a 

a 

CM 

r«. 

CM 

OI 

CD 

C3 

T- 

a 

CO 

a 

EO 

CM 

a 

o 

in 

o 

01 

CM 

rx 

a 

CO 

a 

a 

yf 

a 

y^ 

cvi 

CD 

CO 

CO 

r*. 

O 

a 

a 

o 

a 

a 

a 

a 

a 

a 

EO 

CD 

CO 

rx 

lx 

EO 

a 

o 

a 

a 

a 

a 

CO 

rv 

fx. 

rx 

rx 

rx 

fx 

a 

a 

a 

a 

fx 

fx 

rx 

rx 

o 

CM 

■r. 

o 

y- 

CM 

a 

a 

'M- 

a 

a 

a 

y~ 

a 

CM 

09 

O 

CM 

■M* 

o 

O 

■»- 

a 

a 

a 

rx 

a 

CO 

CD 

CM 

in 

rx. 

CM 

CM 

rx 

fx 

rx 

a 

o 

a 

a 

a 

CD 

in 

OI 

CO 

CO 

CO 

a 

o 

a 

a 

a 

a 

a 

CD 

o 

o 

CM 

CD 

rx 

to 

fx 

o 

o 

a 

a 

a 

CO 

K 

rx- 

rx 

rx 

rx 

rx 

rx 

a 

a 

a 

a 

rx 

rx 

rx 

|X 

in 

•M- 

CM 

rx 

a 

rx 

CM 

'M* 

CM 

a 

o 

a 

CM 

rx 

CO 

r- 

CO 

o 

CD 

CM 

a 

a 

a 

V 

a 

a 

a 

a 

OI 

■M- 

in 

CM 

a 

a 

fx 

CO 

a 

CO 

fx 

a 

rx 

ed 

CD 

CM 

o 

a 

-M- 

a 

'M- 

O 

fx 

a 

a 

o 

a 

a 

a 

r- 

CD 

o 

in 

<o 

a 

a 

O 

O 

o 

o 

a 

a 

a 

a 

co 

CO 

CO 

jhx 

rx 

rx 

irx 

fx 

a 

a 

a 

a 

rx 

rx 

fx 

fx 

Tm 

"T^ 

To 

"a 

_ 

“ 

a 

To 

"a 

To 

"d 

CD 

' 

CD 

CO 

■vr 

a 

Cl 

OI 

OI 

>  CO 

o 

o 

m; 

d 

ed 

cvi 

EO* 

oi 

V- 

y^ 

oi 

y^ 

a 

y^ 

CO 

a 

d 

cvi 

in 

CO 

CD 

CM 

rx 

a 

a 

o 

a 

a 

•M- 

o 

a 

CM 

CO 

CO 

O 

CO 

CO 

in 

a 

a 

a 

a 

a 

a 

a 

a 

CO 

CO 

CO 

rx. 

rx 

rx 

h. 

rx 

rx 

fx 

rx 

rx 

fx 

fx 

rx 

rx 

"d 

Tn 

"d 

To 

a 

"oj 

"a 

"d 

"a 

“ 

"d 

CO 

rx. 

rx. 

in 

Cvj 

a 

a 

a 

rx 

O) 

a 

cd 

CM 

y^ 

cvi 

cvi 

oi 

d 

V 

d 

y^ 

cvj 

r>i 

cd 

CD 

to 

CO 

CO 

CD 

rx. 

'M' 

a 

a 

CM 

a 

a 

a 

a 

o 

CD 

r«- 

CO 

O) 

CO 

CO 

a 

a 

a 

O 

a 

a 

a 

a 

CO 

co 

CO 

CO 

rx 

rx 

rx 

rx 

rx 

rx 

a 

fx 

IX 

fx 

fx 

fx 

"d 

To 

"d 

"a 

To 

"a 

Tn 

Tn 

"d 

CO 

CD 

CO 

CO 

CD 

OI 

a 

CO 

a 

CO 

T- 

'*■ 

a 

o 

V* 

OI 

V 

d 

a 

oi 

a 

a 

yJi 

a 

d 

a 

a 

oi 

CO 

CO 

CD 

o 

rx 

a 

rx 

rx 

fx 

'M- 

a 

■M- 

CO 

rv 

CD 

CD 

eo 

« 

a 

't 

a 

a 

o 

a 

O 

o 

O 

o 

CO 

CO 

CO 

CO 

rx 

rx 

rx 

rx 

fx 

rx 

a 

rx 

a 

a 

a 

a 

"d 

"d 

To 

CD 

a^ 

CJ 

Tn 

"d 

GO 

“d 

d 

Tm 

"d 

in 

CM 

rx. 

OI 

in 

CO 

a 

CD 

rx 

CO 

a 

CO 

a 

a 

■M- 

a 

ED 

y^ 

cd 

cvi 

d 

V 

d 

CD 

cvi 

y^ 

r*.^ 

cd 

cd 

d 

oi 

CM 

o 

m 

OI 

m 

■M- 

a 

a 

CM 

rx 

a 

a 

rv 

eo 

fXx 

CD 

eo 

•Vf 

a 

a 

a 

o 

a 

o 

o 

o 

CO 

CD 

CO 

CO 

fx 

rx 

rx 

fx 

r. 

a 

fx 

a 

o 

a 

a 

f 

d 

? 

UJ 

UJ 

d 

f 

d 

UJ 

UJ 

d 

5 

o 

? 

UJ 

UJ 

d 

$ 

d 

? 

UJ 

UJ 

« 

o 

OI 

o 

d 

eo 

CD 

o 

d 

a 

a 

o 

d 

a 

a 

o 

d 

a 

o 

o 

O 

o 

a 

o 

o 

a 

o 

o 

a 

f 

p 

o 

o 

UJ 

? 

f 

o 

5 

UJ 

5 

o 

o 

UJ 

5 

$ 

o 

o 

UJ 

d 

m 

d 

OI 

d 

d 

CD 

o 

to 

d 

Oi 

d 

d 

a 

d 

a 

3 

d 

d 

a 

d 

a 

d 

a 

d 

d 

a 

CO 

CO 

(0 

CO 

un 

d 

in 

d 

eo 

■M- 

EO 

o 

o 

O 

O 

« 

CO 

<n 

■ 

3 

I 

1 

■1 

d 

CO 

3 

1 

1 

1 

1 

1 

UJ  I 

ffil 

i  ^ 
§i 

ii 


s 


8 


Coefficients  of  Height  &  Associated  Refractivltes  (Means)  (Meters/N  Units) 


Performance  Comparison  of  Tropospheric  Propagation  Models 


55 


UJ 

QC 

lU 


X 

z 

cc 

lU 

f 

oc 

o 

z 


> 

0 

2 

8 

g 

0 

3 

< 

1 

► 

< 

z 

S 

oc 

< 

z 

CV 

If 

_ 

« 

309.43 

6032.9000 

323.83 

8210.2000 

316.54 

8063.6550 

Ok  a 

0  « 
d  P 

CM  *- 
0  C 
« 

316.77 

7966.2417 

327.34 

8214.5765 

317,14 

8102.5625 

320.89 

7991.5911 

322.90 

8191.0321 

317.20 

7649.2875 

336.55 

7SliP  .'ll 25 

C 

«c 

c 

u 

0 

« 

cc 

c 

N 

360.58  1 

C 

« 

N 

U 

« 

o 

o 

o' 

0 

7911.9667 

328.55 

8231.9556 

316.10 

8062.8300 

P 

0 

0  C 

^  p 

d  p 

CM  If 

0  0 

p 

0 

0 

d 

0 

0 

0 

0 

0 

d 

0 

0 

p 

326.81 

8054.5703 

324.08 

8048.2900 

322.71 

7772.6938 

332.20 

0 

0 

CM 

0 

0 

0 

CM 

V  0 

■9  0 

M  0 

9  0 

lA 

CM 

0  ^ 

N  C 

W  « 

PD 

N 

h- 

« 

p 

0 

p 

PJ 

p 

0 

0 

0 

« 

0 

0 

d 

p 

0 

c 

a 

o 

<0 

« 

0 

o 

« 

7585.5833 

341.28 

8162.6222 

r-  0 
0  ^ 
»  fsl 

ry  o 
P* 

328.19 

7895  4000 

0 

p 

CM 

CM 

0 

p 

0 

0 

0 

0 

0 

0 

0 

p 

CM 

F-  CM 

P 

»  CM 
■M  0 

19  0 

P 

0 

0 

0 

CM 

0 
»  P 

[9  0 

P 

3<'U.3U 

7544.9401 

341.53 

0 

0 

0 

d 

0 

p 

p 

0 

0 

A  0 

r-  0 

9  pj 

9  0 

9  0 

0 

0 

0 

■M  0 

D  0 

3  0 
0 
9  0 
P 

0 

p 

PJ 

« 

0 

0 

0 

0 

0 

s 

p 

0 

0 

0 

p 

0 

lA 

(M 

O 

a 

(D 

« 

o 

«> 

cii 

n 

7297.0500 

352.81 

7837.5000 

O 

O 

A  -M- 
3  0 

"A  CM 
•>  0 
•>  P^ 
N 

333.52 

7776.0022 

CM 

d 

0 

p 

0 

0 

0 

0 

0 

0 

p 

0 

0 

■«.  0 

A  0 

A  0' 

N  »- 
■A  CM 

P 

0 

3k  0 

p  0 

3  CM 

0  0  C 
9  0  < 
P 

0 

M 

1-^0 
B  0  0 

b  d 

9  CM  0 
9  0 

0 

0 

*“ 

0 

d 

0  1 
p 

0 

0 

M  P 

3  or 

r  P 
r  0 

9  0 

0 

0 
0 
B  0 

li  CM 
D  O' 
9  0 
0 

0 

P 

CM 

0 

0 

~C 

P 

0 

0 

0 

380.94  1  387.26  1 

«D 

lA 

K 

O 

«> 

e» 

CM 

O 

« 

1.2999 

357.37 

7325.5222 

K 

O  0 

0  ^ 

D  r^ 

n  CM 

9  0 

K 

0 

0 

«-  m 
n  W 
v  0 
0 

P'^ 

0 

V 

0 

0 

p 

0 

p^ 

M  1 
P 

P 

0 

0 

■  •“ 

M  *-  < 

0  c9  r 
9  0 

9  *-  C 
P 

0 

9- 

»  ^  k 

M  0  C 

PA  C 
D  0  U 
9  C 

P 

CM 

0 

B  0  F- 
M  0  0 

A  0 
f>  0  P 
9  0  0 

0 

CM 

0 

0  1 

0'  c 

CM  F 
0  P 
P 

0 

P 

-  0  C 

n  0  < 

li  0  F 

r  0  k 

9  0  c 

0 

0 

9  0 
3  kO 

-  kb 

0  F- 
9  P 

0 

C 

0 

0 

0 

Ta 

CM 

0 

0 

d 

CM 

0 

0 

« 

o 

« 

p> 

oi 

o 

aa 

«-  1 

0 

d 

CM 

n 

3  0 

3  0 

A  V  0  < 

A  o  o  r 
d  0  ^ 
A  0  *-  u 
0  CM  c 

h» 

0 

A  0 

-  0 
r  d 

A  0 

9  ^ 

P* 

0 
0 
A  »- 
9  •- 
D  ^ 
r  CM 
•>  P- 

0 

PA 

0 

0 

0 

0 

0 

0  < 

0  t 
0  C 
0  C 

P 

CM 

r  0  p 
A  0  c 

B  PA  S 
9  ^  k 
9  C 

p 

0 

P 

>»  P  <1 

M  0  « 

B  kB  S 

B  0  1 
9  P 

P 

0 

0 

t  F-  0 

A  CM  Y 

B  0'  CM 
r  P  P 

9  0  0 

0 

0 

0  h 
0  C 
PA  0 
0  M 

K  ' 

0 

P 

«  0  k 

M  0  C 

k  0  F 

r  0  e 
9  0  c 

P 

0 
0 
0  CM 

3  F- 

“  PA 
0  CM 
9  P 

0 

CM 

P 

CM 

0 

0 

"0 

0 

p 

0 

0 

0 

0 

0 

60 

«> 

<0 

o 

«A 

(A 

• 

lA 

O 

0 

o 

0 

0 

•f 

n 

>  ^ 

A  0  T-  p 
D  CM 

n  d  0 

0  ^  ^  ^ 
9  0  0  P 
’•  K 

o 

0 

0 

CM 

0' 

0 

0 

pi. 

0 
0 
M  CM 
r  ^ 

-  CM 

r  0 
9  0 

p 

0 

0 

d 

CM 

0 

0 

0 

0  « 
CM  P 

0  C 

0  P 
0  P 
P 

0 

0 

e  0  p 

-  p^  s 
A  P^  k 
9  0  U 
9  0  P 
P 

0 

••0  0 
»  P  0 

f  d  M 

a  0  M 

9  9-  p 
p 

0 

0 

A  0  0 
k  0  0 
r  0  0 

r  ^  0 

A  *“  0 
P 

0 

0 

0  c 
0  0 

c 

0  M 

p  p 
p 

0 

P 

A  0  S 
k  CM  P 
!  P 

t  0  p 
P 

0 
P 
B  0 
A  0 

F  0 

>  0 
A  0 
0 

0 

d 

0 

0 

a 

0 

0 

0 

0 

0 

0 

0 

CM 

0 

0 

0» 

n 

« 

0 

K 

N 

0 

o 

« 

0  « 

>*  0 

D  0 

A  O  0  U 
A  0  CM  ^ 

^  d  d  u 

A  0  ®  r 

B  0  O  p 

«• 

0 

0 

0 

0 

P^ 

P' 

p 
0 
-  0 
r>  0 
M  d 
9  0 
9  a 
p 

0 

CM 

CM 

0 

0 

0 

0  C 

P  » 

0  f 
P 

CA 

0 

M  CM  U 

■  T  * 

>  0  C 
P  '»  M 

>  0  P 
P 

0 

A  0  h 
k  0  p 
ti  V  u 

r  CM  p 
A  P  p 

7482.5953 

344.06 

0 

0 

0  « 

CM  V 

p  a 
0  0 

S  *- 

0 

0 

A  0  M 

.  10  « 
i  0  F 
1  0  h 
A  0  P 
0 

0 

0 

CM 

0 

0 

0 

0 

0 

kb 

0 

0 

0 

0 

0 

0 

0 

0 

CM 

P 

0 

0 

0 

o 

0 

0 

o 

0 

c 

c 

0  U 

v  U 

^  u 

M 

0  ■ 

A  0 

A  K 

A  O  C 

A  ^  C>  « 
A  d  PA  C 
r  CM  0  A 

A  0  o  r 

0 

0 

0 

>  N  5 

Ok  s 

d  c 
a  c 
Ok  p 
K 

0 
0 
r  0 
D  0 
li  d 
M  0 
A  0 
p 

0 

P. 

0 

0 

P 

CM  « 

PA  a 

0  P 
0  P 
P 

0 

CM 

0  M 

>00 

CM  « 
0  P 
0  P 
P 

0 

0  c 

1  0  P> 
d  M 
1  P  Cl 
0  P 

7700.8807 

330.00 

0 

CM 

0  kC 

P  c 

p  -p 
0  Cl 

0  p 

0 

0 

0 

CM  0 
0  C 

cii  a 
0  If 

^  p 
p 

0 

0 

0 

0 

0 

0 

P 

p 

0 

PA 

0 

0 

"0 

o 

0 

kO 

cii 

'  0 

0 

p 

0 

d 

0 

0 

"w 

0 

0 

n 

cJ 

0 

o 

0 

h 

C 

0  »■ 

0  Cl 

d  r- 

0  c 

323.26 

8257.1111 

313  86 

p«. 

0 

0  a 

p^  p 

d 

n  «■ 

p 

0 

CM 

CM 

0 

0 

0 

0 

0 

P 

0 

0 

0 

P 

0 

»-  cv 

T  ® 
p  cv 

CM  P! 

0  p: 

p 

0 

0  M 

CM  If 

PA  C 

0  r- 
F-  p 
0 

0 

P 

F~  « 
01 

a 

0  «- 
—  0 

0 

8003.8375 

323.01 

0 

«  ^ 

CM  n 
PA  lb 
0  p- 

CM  0 

0 

0 

p 

■  2  * 

t  0  0 

0  p: 
0  0 

^  p: 
p 

■"o 

0 

0 

p 

cii 

0 

p 

p 

0 

0 

0 

0 

"0 

0 

CM 

0 

P 

P 

0 

kb 

0 

0 

FEB 

"S' 

0 

« 

0 

PI 

« 

o  « 

T  ^ 

d  « 

c 

0  C 

« 

321.43 

8175.3656 

316.56 

0 

0 

)  0  M 
>  0  If 

'  0  « 

CM  ^ 

'  »-  r. 

0 

P 

0 

0 

CM 

d 

0 

0 

0 

<D 

ID 

0 

8014.0750 

324.54 

0 

■  0  P3 

! 

’  0  ^ 

1  0 
r  CM  0 

0 

p 

0 

»  0  0 
!  ^  *■ 
'  PA  0 

1  w  0 
0 

0 

0 

P  0 

CM  0 

0  d 

0  CM 

0  0 

0 

0 

P 

0  CM 
0  Ok 
PA  ^ 
CM  F- 
«  0 
0 

0 

0 

CM  0 

0  0 

F-^  0 

p  0 
P  0 

p 

0 

0 

d 

0 

p 

p 

0 

P 

0 

0 

0 

"0 

0 

CM 

0 

CM 

P 

0 

0 

0 

JAN 

'o' 

o 

i 

o 

o 

*-  o 

d  isi 

o  a 

0  o 

0 

321.05 

8188.3111 

317.23 

0  « 

CM  0 

d  p> 

f-  P3 

0 

0 

0 

P 

0 

d 

i 

p 

0 

0 

8036.7250 

324.08 

8275.5323 

313.04 

8147.4600 

319.56 

0 

0 

0  CM 

0  0 

0  d 

0  F- 

0  0 

0 

P 

0 

0  0 
0  0 
d  kb 

0  F- 

*-  0 

0 

7805.6500 

329.63 

"0  ’ 
p 

0 

CM 

P 

P 

P 

0 

0 

0 

lo 

0 

0 

PA 

0 

CM 

P 

0 

0 

cii 

0 

0 

c 

E 

• 

Su 

II  1 

oz  » 

ill 

z  S  z 

ll 

II 

fl 

ti 

II 

!i 

ll 

H 

1 

• 

0 

n 

1 

z 

v 

.C 

J 

z 

• 

• 

n 

3 

a 

C 

o 

t  ^ 

s  8 

o  o 

?  ? 

o  o 

lU 

o 

0 

o 

Ul 

8 

? 

0 

0 

0 

s 

1 

? 

0 

Ui 

s 

0 

UJ 

d 

Ul 

d 

0 

0 

Ul 

8 

? 

0 

0 

5 

d 

0 

? 

8 

0 

? 

0 

Ul 

1 

0 

UJ 

d 

Ul 

d 

0 

0 

Ul 

8 

? 

0 

0 

£ 

d 

0 

? 

8 

? 

0 

• 

• 

u 

3 

« 

z  z 

s  s 

2  2 

z  z 

>  0 

A  V 

z 

S 

o 

Z 

0 

z 

s 

0 

z 

0 

z 

d 

0 

0 

z 

0 

z 

0 

0 

z 

d 

0 

z 

0 

0 

z 

d 

0 

z 

0 

0 

z 

d 

0 

z 

0 

0 

z 

d 

0 

z 

d 

0 

0 

z 

0 

z 

k 

0 

z 

0 

z 

d 

0 

0 

Z 

0 

z 

k 

0 

z 

A 

s 


I 

a 


< 

oc 

o 

< 

o 


1. 

a 

0 

0 

!  0 
p 

0 

c 

If 

0  Cl 
0  Cl 

B  ac 

9  kC 
kC 

C 

a 

e 

C\ 

D. 

tr 

r- 

a 

« 

386.24  1 

3 

C 

a 

0 

0 

CM 

[  PA 

0 

If 

p 

0  p 

lA  CV 
0  0 
[9  Ml 
« 

p 

a 

0 

a 

p 

flO 

0 

0 

0 

« 

« 

0 

0 

0 

0 

CM 

CM 

'  F.; 

p 

0 

0 

M 

■-  « 
•M  0 

B  M 
9  0 

0 

IT 

« 

a 

r 

0 

a> 

p 

0 

0 

0 

0 

0 

0 

g 

0 

0 

0 

0 

0 

0 

p 

0 

O 
0 
9  CM 
B  P 
b  kE 

9  0 

0 

« 

a 

a 

<9 

0 

0 

0 

0 

p 

0 

p 

CM 

Cii 

0 

0 

AUG 

CM 

CM 

CM 

at  c 
P 

0 

0 
0 
9  P 

3  F- 

••  p 
-  0 
9  0 

0 

0 

0 

0 

a 

0 

0 

0 

0 

kb 

0 

p 

0 

d 

p 

0 

JUL 

0 

0 

0  < 

0  M 

d  a 

0  h 

P  c 

0 

0 
0 
3  0 

tA  0 
b  MA 
-  P 
9  0 

0 

0 

0 

0 

0 

0 

0 

0 

kO 

p 

0 

p 

0 

CM 

0 

d 

p 

0 

Z 

3 

0 

0 

0  C 

F-  P 

P  c 

0  a 

0  p 
0 

0 

CM 

at 

A  0 
0 

0 

0 

d 

■0 

0 

0 

0 

kb 

0 

0 

0 

p 

d 

0 

0 

MAY 

0 

0 

0  « 

F-  P 

0  F 
0  0 
0  P 

0 

0 

0 

PA 

0 

0 

0 

0 

d 

0 

P 

0 

0 

CM 

ab 

0 

0 

0 

d 

0 

0 

APR 

0 

CM 

0  0 

0  ai 

0  Cl 

<M  a 

0  p; 

0 

"  "o 
0 
0 
« 
d 

PA 

0 

0: 

CM 

d 

”0 

kb 

p 

0 

p 

p 

d 

0 

0 

MAR 

P 

0 

P  « 

0  F- 

0  « 

P  PI 

0 

'  "o 
0 
CM 
0 

0 

« 

0 

0 

0 

0 

d 

a 

0 

"cm 

CM 

CM 

0 

P 

0 

0 

0 

p 

0 

0 

a 

Bf 

at 

0  0 

at  0 
ab  PA 

0  P 
0  0 
0 

' 

P 

I  0 
t  P 

i  0 

)  0 
0 

d 

0 

0 

"0 

0 

CM 

0 

d 

0 

p 

0 

0 

0 

0 

0 

0 

JAN 

6862.5225 

373.04 

.^. 

0 

0 

0 

0 

0 

0 

0 

p' 

0 

0 

"d ' 

CM 

CM 

0 

P 

0 

0 

0 

0 

c 

^  E 

? 

1 

¥ 

"e 

w 

U 

xz 

Z 

• 

•C 

z 

• 

0 

3 

0 

c 

0 

_l 

Ul 

8 

0 

UJ 

d 

UJ 

d 

0 

0 

UJ 

8 

? 

8 

0 

? 

0 

ki 

• 

0 

3 

• 

..J 

z 

0 

m 

Sa 

z 

0 

0 

m 

0 

z 

MA 

0 

n 

MA 

o 


X 

o. 

</> 

s 


z 

oc 

UJ 


3 

o 

(/> 


f  E 

O  z 


— i - 355,79  I  351.80  I  348.80  j  345.40  I  342,95  I  342.06  I  342.64  I  34^.68  I  347.28  I  349?^ 


Appendix  B 

AVERAGE  ANNUAL  SURFACE  REFRACTIVITIES  AND  MONTHLY 
DEVIATIONS  FOR  10  AREAS  OF  INTEREST 


Three  databases — ^European  Center  for  Medium-Range  Weather  Forcast  (ECM),  High-Resolution 
Analysis  System  (HIRAS),  and  Medium-Range  Frequency  (MRF) — have  been  compared  hourly  with 
mean  and  standard  deviation  for  monthly  deviations  from  average  annual  refractivity  for  10  areas  of 
interest. 
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Appendix  C 

REFRACTIVITIES  FOR  HIGH  ALTITUDE  AREAS 
FOR  29  AREAS  OF  INTEREST  BY  SEASONS 

Refractivities  for  high  altitude  areas  above  1000  m  from  the  mean  sea  level  (MSL)  are  compared  by 
seasons  and  continents  with  tropical  separations. 
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Refractivity  Averages 

Areas-of-Interest  with  High  Altitude  Surface  Pressures 

(Altitude  Sourc«:  MRF  Data;  Rafractlvlty  Sourca;  ECM  Data) 
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Appendix  D 

FIRST  AND  SECOND  ORDER  REFRACTIVITY  GRADIENTS 
FOR  12  AREAS  OF  INTEREST 


First-  and  second-order  refractivity  gradients  are  compared  for  12  areas  of  interest  by  months  with 
European  Center  for  Medium-Range  Weather  Forecast  (ECM)  and  High-Resolution  Analysis  System 
(HIRAS)  databases.  It  is  noticed  that  the  second-order  refractivity  gradients  are  good  sources  for  distinction 
from  normal  refractivity,  and  they  are  easily  comparable  with  other  areas,  months,  and  seasons. 
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Junho  Choi 


DATABASE:  ECM  First  Order  Refractivity  Gradients  (lOOOmb  -  850mb) 


A01 

JAN 

FS 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SEP 

OCT 

NOV 

DEC 

CAN 

47.82 

47.83 

48.26 

50.36 

53.47 

56.86 

60.63 

60.82 

56.62 

52.12 

48.14 

47.72 

CAM 

83.98 

83.88 

85.63 

87.77 

88.12 

89.49 

92.37 

92.27 

89.89 

86.99 

86.41 

85.50 

AMFOR 

81.25 

80.47 

82.11 

83.84 

82.74 

81.63 

78.83 

78.80 

79.65 

81.64 

81.69 

81.98 

SAF 

72.63 

73.31 

74.01 

74.63 

72.11 

67.96 

66.05 

65.43 

65.77 

67.76 

69.35 

71.82 

SAH 

46.36! 

44.81 

45.44 

49.78 

56.74 

61.00 

64,23 

66.92 

65.16 

58.73 

54.10 

49.58 

AUS 

69.67 

71.02 

69.94 

69.49 

67.96 

66.93 

64.50 

63.34 

64.26 

65.82 

67,65 

69.82 

SEAS1 

83.60 

84.33 

86.94 

93.00 

94.07 

91.37 

89.34 

89.03 

89.84 

89.42 

86.28 

83.84 

SEAS2 

84.03 

85.52 

89.42 

92.07 

91.90 

89.59 

88.55 

88.55 

88.81 

88,81 

87.14 

84.43 

GOBI 

43.51 

43.48 

44.83 

47.32 

51.59 

55.43 

56.85 

56.38 

51.58 

47.42 

45.74 

44.27 

47.70 

48.01 

49.32 

53.67 

57.13 

60.68 

60.56 

57.87 

55.791 

51.68 

49.34 

48.57 

SB 

50.45 

50.05 

48.91 

47.60' 

47.78 

52.07 

54.66 

53.95 

49.84 

46.12 

47.25 

49.17 

NAK 

46.71 

48.00 

49.29 

51.30 

51.11 

50.59 

48.32 

52.82' 

57.70 

55.55 

51.95 

48.33 

DATABASE:  ECM  First  Order  Refractivity  Gradients  (850mb  -  700mb)) 


AOI 

JAN 

FS 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SS> 

OCT 

NOV 

DEC 

CAN 

44.66 

44.34 

43.90 

44.37 

47.69 

51.89 

55.88 

55.19 

50.97 

47.67 

45.69 

44,87 

CAM 

67.43 

68.03 

67.66 

67.51 

71.67 

70.89 

68.95 

68.78 

69.67 

70.13 

71.06 

68.89 

71.07 

71.66 

71.10 

71.14 

72.35 

70.24 

68.48 

68.84 

69.45 

70.81 

72.05 

71.79 

SAF 

59.06 

60.32 

60.29 

58.90 

55.87 

54.41 

52.80 

52.32 

53.95 

54.36 

56.34 

58.09 

SAH 

48.50 

45.34 

44.05 

46.04 

48.62 

52.19 

57.32 

60.51 

56,60 

53.34 

51.13 

49.28 

AUS 

58.54 

59.64 

59.16 

56.65 

56.94 

53.50 

52.49 

52.61 

51.21 

52.37 

54.55 

56.22 

SEAS1 

69.26 

68.50 

68.71 

68.26 

69.21 

69.60 

70.63 

69.82 

69.13 

70.15 

69.85 

70.11 

70.09 

68.94 

67.63 

69.41 

70,96 

70.75 

68.80 

67.69 

67.32 

68.76 

70,52 

69.92 

44.80 

44.36 

44.55 

45.46 

47.99 

50.03 

52.47 

52.76 

49.02 

46.08 

45,26 

45.50 

44.01 

44.32 

45.91 

49.23 

52.42 

54,68 

53.42 

49.14 

47.16 

45.65 

45.35 

SB 

43.97 

43.71 

42.97 

42.59 

43.83 

47.31 

51.13 

49.95 

46.50 

44.09 

44.00 

43.98 

NAK 

45.53 

45.51 

45.57 

45.58 

47.80 

50.70 

55.87 

56.74 

51.51 

48.93 

47.49 

46.49 

DATABASE:  ECM  Second  Order  Refractivity  Gradients  [(lOOOmb  -  850mb)  -  (850mb  -  700mb)] 


AOi 

JAN 

FS 

MAR 

APR 

MAY 

JUN 

JUL 

AUG 

SB> 

OCT 

NOV 

DEC 

CAN 

3.16 

3.48 

4.36 

5.99 

5.78 

4.97 

4.75 

5.63 

5.65 

4.45 

2.45 

2,85 

CAM 

16.55 

15.84 

17.96 

20.26 

16.45 

18.60 

23.42 

23.50 

20.22 

16.86 

15.35 

16,62 

AMFOR 

10.18 

8.81 

11.00 

12.70 

10.38 

11.39 

10.35 

9.97 

10.20 

10.83 

9.64 

10.19 

SAF 

13.57 

12.99 

13.73 

15.73 

16.23 

13.55 

13.25 

13.10 

11.82 

13.40 

13.01 

13.73 

SAH 

-2.14 

1.40 

3.74 

8.12 

8.81 

6.91 

6.41 

8.57 

5.40 

2.97 

0.30 

AUS 

11.13 

11.38 

10.78 

12. as 

11.01 

13.43 

12.01 

10.74 

13.05 

13.44 

13.10 

13.61 

SEAS1 

14.34 

15.83 

18.23 

24.74 

24.86 

21.77 

18.71 

19.21 

20.71 

19.27 

16.43 

13.73 

SEAS2 

13.94 

16.57 

21.79 

22.65 

20.93 

18.84 

19.75 

20.86 

21.49 

20.05 

16.62 

14.51 

-1.30 

-0.88 

0.28 

1.86 

3.60 

5.40 

4,38 

3.62 

2.56 

1.35 

0.49 

-1.24 

3.22 

4.01 

5.00 

7.76 

7.90 

8.26 

5.87 

4.45 

6.65 

4.52 

3.69 

3.22 

SB 

6.48 

6.34 

5.95 

5.01 

3.95 

4.77 

3.53 

4.00 

3.34 

2.03 

3.25 

5.19 

NAK 

1.18 

2.49 

3.72 

5.71 

3.31 

-7.55 

-3.92 

6.19 

6.62 

4.46 

1.84 
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DATABASE:  HIRAS  First  Order  Refractivity  Gradients  (lOOOmb  -  850mb)  6.50 

API  I  I  JAN  I  Fg  I  MAR  I  APR  I  MAY  I  JUN  I  JUL  I  AUG  I  SB»  I  OCT  I  NOV  I  DEC 

CAN _ 00  46.54  45.57  44.65  45.38  47.28  51.37  54.66  55.84  53.90  49.86  47.37  46.77 

CAN _ 06  46.39  46.08  46.10  47.21  49.98  54.58  58.99  59.52  55.92  51.09  47.75  46.66 

CAN _ 12  47.25  47.38  47.56  48.74  51.90  55.62  59.40  60.74  56.72  51.58  47.40  47.16 

CAN _ 18  46.30  46.09  46.14  46.68  48.90  52.51  56.67  57.53  54.43  50.25  47.40  46.61 

CAM _ 00  82.62  82.81  84.83  63.87  81.81  82.00  81.98  82.31  82.08  81.94  81.70  82.37 

CAM  06  83.01  84.13  85.40  85.71  86.04  87.19  88.04  87,04  85.95  84.04  83.17  R3  Afl 

CAM _ 12  80.13  60.63  83.13  83.23  83.14  82.88  83.60  83.12  81.50  80.61  75.10  80.10 

CAM _ 18  74.44  75.52  76.28  76.36  77,30  78.96  79.45  78.70  78.04  75.49  75.10  74.75 

AMFOR  00  68.12  67,79  69.79  70,56  73.79  76.58  77,94  79.40  77.99  75.75  73,55  69.96 

AMFOR  06  66.28  65.83  66.40  67.13  69.98  74.13  76.92  78.87  77.42  75.33  72.15  68.72 

AMFOR  12  70.83  71,14  72.76  72.08  73.15  74.96  75.33  76.47  76.15  75.60  62.08  71  Qfl 

AMFOR  18  58.82  58.35  58.89  59.41  61.54  65.04  65.60  66.08  64.72  63.50  62.03  59,79 

SAF _ 00  73.94  71,66  72.63  69.61  67.09  61.74  62.56  63.16  65.88  65.38  69.19  70,60 

SAF _ 06  67.37  65.46  66.83  64.01  63.33  59.20  60.95  61.78  62.77  62.04  64  48  64  SB 

SAF _ 12  72.04  71.03  71,45  69.06  67.60  60.39  61.76  62.56  65.98  65.54  62.90  69.70 

SAF _ 18  62.95  60.92  62.73  61.03  60.83  57,72  58.74  59.73  61.00  60.64  62.90  61.39 

SAH _ 00  71.43  71.02  69.80  69.36  74.65  80.65  82.83  83.49  65.56  77.72  71.97  69.88 

SAH _ 06  72.99  71,71  70.11  67.88  72.32  77.82  60.67  81.97  83.24  76  93  71  an  71  flP 

SAH _ 12  67.50  65.95  64.47  62.30  66.37  71.93  72.16  74.03  75.66  69.97  63  48  66  ng 

SAH _ 18  70.80  70,05  66.30  65.05  68.22  72.11  74.21  76.69  78.40  72.66  69.46  69.28 

AUS _ 00  67.22  67.82  66.20  62.25  59.62  59.01  58.19  58.31  60.67  63.74  65.30  68.12 

AUS _ 06  58.40  58.40  57.22  55.55  55.62  57,14  56.78  56.43  57.57  58.92  59.65  59.38 

AUS _ 12  62.96  63.28  62.51  59.60  59.15  59.78  59.19  59.20  60.18  62.75  59  38  fins? 

AUS _ 18  56.69  57.86  57.89  57.00  56.93  58.44  57.87  58.10  58.43  59.25  59.33  58.70 

SEAS1  00  74.11  76.80  79.50  81.78  80.32  79.58  77.95  76.88  77.13  76.98  74.65  71.86 

SEAS1  06  61.29  63.80  65.87  68.23  67.89  70.11  68.74  68.32  68.27  67.35  64.87  60.72 

SEAS1  12  70.14  72.26  74,48  77.27  77,14  76.55  75.21  75.03  74.66  74.14  69.78  68.86 

SEAS1  18  66.76  69.89  73.05  75.77  75.21  75.88  74.10  73.53  74  07  72.83  fiQ  7H  fiS  Q1 

SEAS2  00  75.02  76.71  77.43  79.26  79.07  79.13  78.25  78.51  78.64  77.42  75.47  74.29 

SEAS2  06  61.43  63.12  63.76  65.73  67.34  69.75  68.27  68.72  67.93  66.34  64.03  61.21 

SEAS2  12  73.33  75.72  77  01  79.40  79.54  79.95  79.48  79.40  78.77  77.52  72.11  73.26 

SEAS2  18  68.43  70.45  72.12  75.16  76.81  78.37  77.89  77.97  77.77  75.50  72.11  68.41 

GOBI _ 00  46.05  44.72  43.91  45.97  51.97  59.37  64.01  65.60  58.77  52.61  46.94  46.13 

GOBI _ 06  44.59  43.39  42.46  43.52  47.77  53.12  59.54  61.45  54.06  50.02  46.25  45.26 

GOBI _ 12  44.28  42.43  42.07  43.90  51.45  59.20  65.79  69.95  62.21  52.98  46.89  44,96 

GOBI _ 18  45.00  43.44  42.97  44.48  48.78  54.52  59.84  62.30  56.47  52.33  46.89  45.63 

EURAS  00  47.25  47,54  48.99  52.01  56.46  63.84  67.02  64.54  59.00  53.72  49.78  48.29 

EIWAS  06  46.74  46.76  47.93  49.43  51.68  58.47  62.04  60.36  55.73  51.66  49.23  47.91 

EURAS  12  46.31  45.90  45.88  45.51  47,28  53.65  56.41  55.51  51.61  48.74  49,13  47.55 

EURAS  18  46.79  46.74  47.64  ,49.09  52.14  58.13  61.87  60.80  55.89  51.71  49,13  47.90 

SB _ 00  50.04  49.09  47.11  45.44  45.37  49.43  54.20  52.73  48.45  44.69  47.05  48.78 

SB _ 06  47.88  47,04  45.19  43.86  44,04  47.19  51.54  50.36  46.97  44.18  45.61  46.88 

SB _ 12  49.66  48.70  46.03  44.13  44,68  48.76  52.71  52.07  48.07  44,44  46.14  48.77 

SB _ 18  48.64  47.64  45.91  44.70  45.84  50.59  55.77  53.92  48.84  44.60  46.14  47.56 

NAK _ 00  45.46  45.86  45.65  47.39  49,56  52.60  54.77  56.92  53.62  50.59  47.45  46.11 

NAK _ 06  44.95  45.16  45.27  47.00  49.33  53.30  57.01  57.79  53.05  49.60  46.76  45.59 

NAK _ 12  45.42  45.60  45.74  47.54  50.20  52.99  54.72  57.57  54,06  50.66  46.81  45.95 

NAK  I  18  44.97  45.041  45.28  47.39  49.84  53.60  56.60  58.28  53.34  49.92  46.81  45.83 
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DATABASE:  HIBAS  First  Order  Refractivity  Gradients  (850mb  >  700mb)) 


AOI 


CAN 


CAN 


CAN 


CAN 


CAM 


CAM 


CAM 


CAM 


AMFOR 


AMFOR 


AMFOR 


AMFOR 


SAF 


SAF 


SAF 


SAF 


SAH 


SAH 


SAH 


SAH 


AUS 


A(^ 


AUS 


AUS 


SEAS1 


SEASI 


SEAS1 


SEASI 


SEAS2 


SEAS2 


SEAS2 


SEAS2 


GOBI 


GOBI 


GOBI 


EURAS 


EURAS 


EURAS 


EURAS 


se 


S6 


SB 


SB 


NAK 


NAK 


K 


NAK 


APR 


43.55 


44.99 


44.60 


43.04 


65.78 


70.07 


69.84 


64.42 


71.56 


80.93 


70.59 


70.70 


55.10 


59.79 


53.87 


61.60 


44.86 


45.01 


41.93 


43.61 


57.01 


61.57 


60.65 


63.68 


72.84 


71.88 


70.30 


79.76 


72.96 


72.63 


75.42 


83.07 


43.46 


40.99 


38.00 


40.55 


46.41 


44.07 


44.29 


MAY 


45.27 


47.91 


46.74 


44.02 


67.18 


72.52 


70.19 


66.35 


71.00 


79.98 


70.28 


69.05 


52.75 


56.76 


51.55 


59.26 


46.10 


46.06 


42.77 


44.40 


57.70 


59.51 


59.08 


60.21 


72.96 


72.76 


71.24 


79.63 


73.37 


73.24 


75,36 


82.53 


45.78 


43.04 


40.01 


42.55 


48.49 


44.68 


45.44 


49.80 


44.44 


43.56 


44.18 


45.09 


48.16 


48.12 


48.14 


47.96 


JUN 


49.20 


52.72 


50.40 


46.79 


67.91 


73.54 


71.17 


67.97 


67.72 


75.48 


67,76 


65.17 


52.18 


54,04 


50.93 


55.74 


45.69 


46.78 


43.56 


44.63 


56.63 


56.81 


56.85 


57.10 


71.76 


71.51 


70.39 


77.36 


71.25 


70.51 


72.83 


79.65 


51.96 


47,49 


44,67 


46.43 


51.97 


47.04 


48.12 


54.44 


46.89 


45.60 


46.73 


48.78 


49.88 


50.67 


50.15 


50.37 


JUL 


52.42 


56.51 


53.39 


49.16 


68.48 


73.90 


71.06 


67.20 


64.40 


71.78 


64.70 


61.20 


51.63 


53.42 


50.44 


55.80 


49.48 


50.18 


46.62 


47.31 


54.93 


54.79 


55.59 


55.54 


69.96 


69.41 


69.11 


75.32 


69.69 


68,67 


71.06 


77.94 


57.77 


54.68 


49.46 


51.51 


54.70 


49.63 


50.71 


57.31 


50.23 


48.36 


49.94 


52.89 


54.13 


55.03 


54.29 


54.64! 


AUG 


52.60 


56.05 


53.18 


48.81 


68.51 


73.49 


70.47 


67.20 


64.59 


72,83 


64,81 


60.45 


52.35 


54.44 


51.52 


56,51 


51.74 


51.67 


48.59 


49.34 


52.89 


53.52 


54.27 


54.95 


70.38 


70,20 


69.15 


75.22 


69.16 


68.07 


70.74 


77.02 


57.06 


54.28 


48.88 


51,34 


54.05 


49.65 


50.04 


56.45 


49.36 


47.41 


49.10 


50.98 


54.92 


55.67 


54.73 


55.06 


S^  I 


49.09 


51.10 


49.92 


46.71 


67.92 


72.43 


69.42 


67.02 


66.07 


75.11 


65.75 


62.30 


53.21 


56.03 


51.88 


57.04 


49.10 


49.50 


45.89 


47.24 


52.84 


54.72 


55.22 


57.10 


70.54 


70.29 


69.47 


75.66 


69.05 


68.22 


70.82 


77,96 


50.72 


47.29 


44.57 


47.28 


50,41 


46.90 


47.55 


51.69 


46.30 


44.99 


46.24 


46.78 


51.83 


51.91 


52.28 


51.53 


cxrr 


46.12 


46.97 


46.63 


45.13 


67,78 


71.27 


69.09 


65.88 


67.72 


77.25 


66.96 


64.88 


55,31 


57.77 


54.26 


58.63 


47.69 


48.98 


44.22 


46.54 


52.30 


55,35 


55.90 


58.54 


71.01 


69.92 


70.05 


75.83 


70.08 


69.85 


72.72 


79.51 


46.64 


43.87 


42.46 


45.22 


47.02 


45.13 


45.20 


47.59 


44.21 


43-80 


44.10 


44.17 


48.75 


48.35 


49.00 


48.331 


NOV  I 


44.85 


45.02 


44.56 


44.56 


68.39 


70.91 


66.18 


66.18 


69.42 


78.55 


67.24 


67.24 


55,73 


59.28 


60.96 


60.96 


49.16 


49.52 


47.29 


47.29 


54.89 


58.75 


61.66 


61.66 


70.48 


68.95 


75.31 


75.31 


71.91 


71.62 


80.48 


80.48 


45.94 


43.57 


44.18 


44.18 


45-34 


44.81 


45.52 


45.52 


44.01 


44.17 


44.03 


44.03 


47.01 


46.34 


46.14 


46.14 


DEC 


44.35 


44.24 


44.33 


43.96 


66.33 


67.85 


69.99 


63.61 


69.95 


78.27 


68-73 


68.74 


58.17 


62.50 


57.70 


63.86 


49.01 


49.09 


46.05 


46.64 


56.01 


60.80 


60.94 


64.11 


71.68 


68.73 


70.51 


74,97 


72.^5 


71,68 


74.86 


79.82 


46.65 


44.49 


45.15 


44.55 


45.22 


44.93 


44.97 


45.20 


44.31 


44.63 


44.20 


44.33 


45.94 


45.41 


45.80 


45.33 
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DATABASE;  HIRAS  Second  Order  Refractlvlty  Gradients  [(lOOOmb  -  850mb)  -  (850mb  -  TOOmb)] 


CAN 


CAN 


CAN 


CAN 


CAM 


CAM 


CAM 


CAM 


AMFOR 


AMFOR 


AMFOR 


AMFOR 


SAF 


SAF 


SAF 


SAF 


SAH 


SAH 


SAH 


SAH 


AUS 


AUS 


AUS 


AUS 


SEAS1 


SEAS1 


SEAS1 


SEAS1 


SEAS2 


SEAS2 


SEAS2 


SEAS2 


GOBI 


GOBI 


GOBI 


GOBI 


EURAS 


EURAS 


EURAS 


EURAS 


SB 


SB 


SB 


SB 


NAK 


NAK 


NAK 


NAK 


Appendix  E 

ELECTRONIC  COUNTERMEASURES  AND  HIGH-RESOLUTION 
ANALYSIS  SYSTEM  DATABASES 
WORLDWIDE  TIME  DELAYS  BY  SEASONS  AND  MODELS 


Global  time  delays  for  electronic  countermeasures  and  high-resolution  analysis  system  databases  are 
compared  for  each  model  by  seasons  and  elevation  angles  from  the  horizon  (0°)  to  the  zenith  (90°).  Note 
that  the  stratified  model  represents  the  empirical  data  supplied  by  the  Environmental  Technical  Applications 
Center  (Scott  Air  Force  Base). 
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0.1  I  359.4677I  282.2441  284.18631 


253.2324 

1  204.4954 

206.0520 

140.9563 

119.3441 

120.3625 

81.4883 

82.2334 

70.4243 

61.2371 

61.7209 

50.6533 

44.2231 

44.6220 

22.9834 

23.2001 

15.7738 

15.9439 

11.7751 

10.4652 

10.4315 

9.0250 

8.1914 

7.9969 

50.0 


90.0 


TIME  DELAY:  ECM  Data  Global  Sprin 


345.7371 


277.33661 


219.9158 


115.6715 


76.4959 


56.5391 


40.6301 


20.9653 


14.3706 


9.4026 


7.2066I 


271.66481 


232.7362 


196.3543 


116.6830 


79.9267 


59.9972 


43.7459 


16.15901 


11.1379 


9.2428 


283.8239 


246.8413 


210.9880 


128.6369 


89.3642 


67.0020 


47.5619 


14.2234 


8.2859 


5.8546! 


363.2622 


311.4307 


259.4160 


145.5854 


97.8481 


72.9238 


52.4780 


18.6635 


12.2065 


9.3558 


361.66621 


307.51521 


141.3482 


94.7857 


70.57561 


26.2728 


20.0 


30.0  18.0397 


50.0  11.7977 


90.0 


TIME  DELAY:  ECM  Data  Global  Summer 


282.8062 


242.8188 


204.3608 


119.1655 


81.3153 


61.0267 


44.1317 


22.8692 


15.7572 


10.3135 


286.1367 


246.0119 


207.2844 


121.0511 


82.7132 


62.0867 


44.8900 


23.3413 


16.0412 


10.4953 


347.7041 


278.9719 


221.0621 


115.8339 


76.5578 


56.6577 


40.6878 


20.9332 


14.3775 


9.3957 


7.1956 


271.66481  285.75261 


196.3543  212.4218 


116.6830  129.5110 


79.9267 


59.9972 


43.7459 


23.2593 


16.1590 


11.1379 


9.2428 


89.97151 


67.45731 


22.6684 


14.32001 


365.0867 


312.7799 


260.3187 


145.8127 


97.9306 


72.9627 


27.1830 


18.6661 


11,92151 


9.1375 


10.4310 


7.9780 


7.20321 


TIME  DELAY:  ECM  Data  Global  Fall 


43.7459 


23.2593 


16.1590 


11.1379 


9.2428 


8.4221 


5.95081 


12.2534 


9.3920 


Performance  Comparison  of  Tropospheric  Propagation  Models 


TIME  DELAY:  HIRAS  Data  Global  Winter  OOOOHrs 


Junho  Choi 


TIME  DELAY:  HIRAS  Data  Global  Sprinq  OOOOHrs 


Hopfield  Goad 

(ns)  (ns) 


0.1  349,7325 


- - 


299.5528 


277.9101 


258.4271 


249.4374 


181.5074 


140.0251 


224.5822 


209.8807 


203,15141 


EBEBBI 


204.4141 


152.4465 


119.5270 


269.6415 


245.2284 


224.8639 


215.4158 


150.8816 


114.6205 


91,4804 


250.8466 

232.7362' 


255.7520 
238.9236  ' 


328.3719 

304.6138 


147.6423 


116.6830 


95.3896 


nHHFni _ 

i^kee^hhisedi^kek^bkees^ 

1 1 H  ii  I  I II  mil 


55.6928 


10.0  50.4668 


26.1334 


TIME  DELAY:  HIRAS  Data  Global  Sprinq  OGOOHrs 


Elav  Ang 


60.9104 


54.06631 


48.5442 


44.0354 


22.8241 


54.4074 


48.8654 


44.3309 


23.04981 


49.7121 


44.4903 


48.0750 


43.74591 


23.2593 


156.0184 


124.5108 


102.5479 


74.3360 


46.0363 


21.7932 


186.2770 


144.1173 


116.4212 


57.6125 


52.2188 


27.0637 


EEEIQ^I 


279.78291 


281.4504 
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TIME  DELAY:  HIRAS  Data  Global  Fall  OOOOHrs 


Appendix  F 

TIME  DELAYS  AND  ANGLE  ERRORS  OF  42  AREAS  OF  INTEREST 
FOR  SEASONS/ANGLES  BY  HOURS 


Time  delays  of  42  areas  of  interest  are  compared  for  each  model  by  seasons  and  elevation  angles  from 
the  horizon  to  10°  above  the  horizon. 
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using  ECM  Data  and  Hopfiald,  Qoad  and  Exponanlial  Modala 
(Modala  uaa  ECM  turtaea  Data) 
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Time  Delay  (ne) 

November 

using  ECM  Data  and  Hopfiatd,  Qoad  and  Exponantial  Modala 
(Modala  uaa  ECM  Surfaea  Data) 


Angle  Error  (deg) 

February 

using  ECM  Data,  Goad  Modal  and  Exponential  Modal 
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Angle  Error  (deg) 
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Appendix  G 

TIME  DELAYS  AND  ANGLE  ERRORS  FOR  DATABAASES 
AND  SEASONS/ANGLES  BY  HOURS 

Time  delays  and  angle  errors  are  compared  for  10  areas  of  interest  based  on  databases  and  elevation 
angles  from  the  horizon  to  10°  above  the  horizon  by  six  hourly  intervals. 
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Tim©  Delay  (n»)  for  Selected  Areas-of-interest 
ECM,  HIRAS,  MRF  Data  lor  February  1995 
(0000,  0600,  1200  and  1800  Hours) 
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Time  Delay  (ns)  for  Selected  Areas-of-Interest 
ECM,  HIRAS  and  MRF  Data  for  May  1995 
(0000,  0600,  1200  and  1800  Hours) 
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Time  Delay  (ns)  for  Selected  Areas-oMnterest 
ECM,  HIRAS  and  MRP  Data  for  August  1995 
(0000,  0600,  1200  and  1800  Hours) 
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Time  Delay  (ns)  for  Selected  Areas-oMnterest 
ECM,  HIRAS  and  MRF  Data  for  November  1995 
(0000,  0600,  1200  and  1600  Hours) 
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Angle  Error  (degrees)  for  Selected  Areas-of-Interest 
ECM,  HIRAS  and  MRF  Data  lor  February  1995 
(0000,  0600,  1200  and  1800  Hours) 


rcimagiBiii - 


6)  Northsm  Australia.  Tsnsmi  Desert  (NAUS 


renee  MounUiee  (PYRNES 


9)  Tehran.  Iren  (TEHRAN 


10)  Xining,  China  (XINING 


3)  Bangkok.  Thailand  (BANG 


4)  Washington.  D.C.  (DC 


5)  Alasfca  (NA 


6)  Northern  Australia.  Tanami  Deeert  (NAUS 


renee  Mountalne  (PYRNES 


8)  Spokane.  Washington  (SPO 


9)  Tehran,  Iran 


10)  Xining,  China  (XINING 


litiimi-im-a 


ie  (AHAGR 


2)  Amazon  Feraet  (AMFOR 


Thailand  (BANG 


6)  Northern  Auetmlia.  Tanami  Desert  (NAUS 


renee  Mounteine  (PYRNES 


9)  Tehran.  Iran  (TEHRAN 


im  Xining,  China  (XINING 


Imim-L-LL 


rAlgeria  (AHAGR 


2)  Amazon  Forest  (AMFOR 


4)  Washington.  D.C.  (DC 


6)  Northern  Australia.  Tanami  Desert  (NAUS 


renee  Mounteine  (PYRNE 


9)  Tehran.  Iran 


10)  Xining,  China  (XINING 


iBlUJI-i-LL 


r.  Algeria  (AHAGR 


2)  Amazon  Forest  (AMFOR 


3)  Bangkok.  Thailand  (BANG 


6)  Northern  Auetralia.  Tanami  Deeert 


renee  Mountains  (PYRNES 


8)  Spokane.  WMhington  (SPO 


9)  Tehran.  Iran 


10)  Xining,  China  (XINING 


Elevation  Angle  »  0° _ 

^^*^AS  I^F  (2/15/95) 

■  OOQO  0600  1200  1800  0000  0600  1200  1800 


0.3107  I  0.3163 


0.3003 


Elevation  Angle  »  1° _ 

MWAS  MRF  (2/15/95) 

0000  0600  1200  1800  0000  0600  1200  1800 


0.2486 


Elevation  Angie  »  3* 

MRF  (2/15/95) 

^000  0600  1200  1800  0000  0600  1200  1800 


I  ■■irri  **^**^™*™^trma 


jBA«iaaBUJJiBBlBBIiaS!]OZagmMaei^Mn.^aK^| 


_ Elevadon  Angle  «  5* 


H'^AS  MRF  (2/15/95) 

0000  0600  1200  1800  0000  0600  1200  1800 


ISjif^ggBMSigliliaBgiESggEiaglggaameeanwanne 


_ Elevation  Angle  «  10“ 


HWAS  MRF  (2/15/95) 

_  0600  1200  1800  0000  0600  1200  1800 


Performance  Comparison  of  Tropospheric  Propagation  Models 


93 


Error  (degrees)  for  Selected  Areas-oMnterest 
ECM,  HIRAS  and  MRF  Data  for  May  1995 
(0000,  0600,  1200  and  1800) 
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Error  (degrees)  for  Selected  Areas-of-Interest 
ECM,  HIRAS  and  MRF  Data  for  August  1995 
(0000,  0600,  1200  and  1800  Hours) 
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Angle  Error  (degrees)  for  Selected  Areas-of -Interest 
ECM,  HIRAS  and  MRF  Data  for  November  15,  1995 
(0000,  0600,  1200  and  1800  Hours) 


Appendix  H 

TIME  DELAYS  AND  ANGLE  ERRORS  FOR  HOURS 
AND  SEASONS/ANGLES  BY  MODELS 

Time  delays  and  angle  errors  are  compared  for  10  areas  of  interest  by  seasons  and  elevation  angles  by 
tropospheric  models 
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Time  Delay  (ns)  for  10  Selected  Areas-of-lnterest 
MRF,  Hopfield,  Goad  and  Exponential  Model, 
at  0600  Hours 
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Time  Delay  (ns)  for  10  Selected  Areas-of-lnterest 
MRF,  Hopfield,  Goad  and  Exponential  Model 
at  1200  Hours 
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Angle  Error  (degrees)  for  10  Selected  Areas-of-Interest 
MRF,  Goad  and  Exponential  Model 
at  0000  Hours 
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Angle  Error  (degrees)  for  10  Selected  Areas-oMnterest 
MRF,  Goad  and  Exponential  Model 
at  0600  Hours 
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Angle  Error  (degrees)  for  10  Selected  Areas-oMnterest 
MRF,  Goad  and  Exponential 
at  1200  Hours 
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Angle  Error  (degrees)  for  10  Selected  Areas-of-Interest 
MRF,  Goad  and  Exponential  Model 
at  1800  Hours 


Appendix  I 

TIME  DELAYS  AND  ANGLE  ERRORS  FOR  SEASON 
AND  MODELS/ANGLES  BY  HOURS 

Time  delays  and  angle  errors  are  compared  for  10  areas  of  interest  with  seasons  and  models  by  hours 
from  the  horizon  to  10°  from  the  horizon. 
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Tim©  Delay  (ns)  for  Selected  Areas-oMnterest 
MRF,  Hopfleld,  Goad  and  Exponential  Model  for  15  February  1995 
(0000,  0600,  1200  and  1800  Hours) 
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Angle  Error  (degrees)  for  Selected  AreasoMnterest 
MRF,  Goad  and  Exponential  Model  for  15  November  1995 
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Appendix  J 

TIME  DELAYS  AND  ANGLE  ERRORS  FOR  SEASONS 
AND  HOURS/ANGLES  BY  MODELS 

Time  delays  and  angle  errors  are  compared  for  10  areas  of  interest  with  seasons  and  hours  by  models 
from  the  horizon  to  10°  elevation  angles. 


117 


118 


Junho  Choi 


Time  Delay  (na)  for  Selected  Areas-of-Interest 
MRF,  Goad  and  Exponential  Model  for  15  February  1995 
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Time  Delay  (ns)  for  Selected  Areas-oMnterest 
MRF,  Goad  and  Exponential  Model  for  15  May  1995 
(0000,  0600,  1200  and  1800  Hours) 
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MRF,  Goad  and  Exponential  Model  for  15  November  1995 
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Angle  Error  (degrees)  for  Selected  Areas>of‘lnterest 
MRF,  Goad  and  Exponential  Model  for  15  May  1995 
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Angie  Error  (degrees)  for  Selected  Areas-oMnterest 
MRF,  Goad  and  Exponential  Model  for  15  August  1995 
(ooao,  0600,  1200  and  1800  Hours) 
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Angle  Error  (degrees)  for  Selected  Areasof-lnterest 
MRF,  Goad  and  Exponential  Model  for  15  November  1995 
(0000,  0600,  1200  and  1600  Hours) 


Appendix  K 

TIME  DELAY  VARIATIONS  ON  AZIMUTH  BY  SEASONS  AND  HOURS 

Time  delays  are  compared  by  seasons  and  hours  for  azimuth  ange  variations  from  0°  to  360°  in  two 
noticeably  sensitive  areas — ^Teheran,  Iran,  and  Ahaggar,  Algeria. 
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Ahaggar,  Algeria,  (ECM  surface  data — 0°  elevation  angle) 


Ahaggar,  Algeria,  0000  h  (HIRAS  surface  data— 0°  elevation  angle) 
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Ahaggar,  Algeria,  1800  h  (HIRAS  surface  data— 0°  elevation  angle) 


Ahaggar,  Algeria,  February  (HIRAS  surface  data — 0°  elevation  angle) 
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Ahaggar,  Algeria,  November  (HIRAS  surface  data— 0°  elevation  angle) 


Tehran,  Iran  (ECM  surface  data — 0°  elevation  angle) 
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Tehran,  Iran,  0000  h  (HIRAS  surface  data — 0°  elevation  angle) 


Tehran,  Iran,  0600  h  (HIRAS  surface  data — 0°  elevation  angle) 
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Tehran,  Iran,  August  (HIRAS  surface  data— 0°  elevation  angle) 


Tehran,  Iran,  November  (HIRAS  surface  data— 0°  elevation  angie) 


